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PREFACE 


This  report  describes  work  performed  at  the  Electrical  Engi- 
neering Department  cf  the  University  of  Tennessee,  Knoxville,  Ten- 
nessee 37916,  over  the  period  16  October  1973  to  31  December  197^. 
The  work  was  performed  for  the  Air  Force  Systems  Command,  Rome  Air 
Develoj®ent  Center,  Griffins  Air  Force  Base,  New  York  13^1  under 
contract  F30602-7U-C-0037. 

The  overall  effort  related  to  the  definition  of  a rain  back- 
scatter  model  and  cancellation  of  radar  rain  clutter  using  a new 
polarization  method.  The  work  was  conceived  and  guided  hy  Dr.  Peyton 
2.  Peebles,  Jr.,  Associate  Professor  of  Electrical  Engineering.  Most 
of  the  work  on  the  rain  cancellor  was  performed  by  Mr.  R.W.  Rice,  a 
Graduate  Research  Assistant.  Dr.  Peebles  did  most  of  the  work  re- 
garding the  rain  model,  while  Mr.  H.  Sakamoto,  a Graduate  Teaching 
Assistant,  performed  several  literature  searches.  The  final  report 
was  prepared  by  Dr.  Peebles  and  Mr.  Rice. 

The  technical  contract  officer  for  the  Air  Force  was  Mr.  Vincent 
C.  Vannicola.  The  authors  express  their  appreciation  to  Mr.  Vannicola 
for  his  support  of  the  research. 


CONTESTS 


1.0  INTRODUCTION 

1.1  Historical  Development  — — — — 

1.2  Description  of  the  Proposed  Method  — 


2.0  BACKGROUND  THEORY  AND  ASSUMPTIONS 5 

2.1  Basic  Assumptions " — 5 

2.2  Coordinate  Definitions  -- 6 

2.3  Scattering  Theory 9 

2.U  The  Reference  System 13 

2.5  The  Ordinary  System  — 17 

3.0  THEORETICAL  RAIN  MODEL  — 25 

3.1  Drop  Shape " — . 25 

3.2  Drop-Size  Distribution — — 32 

Measured  Data — 33 

Discussion  of  Measurement  Data 38 

Drop-Size  Distribution  Models 1*0 

Summary  of  Drop-Size  Distribution  Models  51 

3.3  Probability  Density  Function  for  Drop-Size  51 

U.Q  ORDINARY  SYSTEM  PERFORMANCE  53 

l+.l  System  Model  -L— — — — — 53 

U.2  System  Clutter  Powers  — 5** 

Ordinary  System — 51* 

Reference  System — — 55 

!*.  3 Clutter  Cancellation  Ratio 55 

Ideal  Rain  Case — 56 

Ideal  System  Case 56 

U.U  Ideal  System  Numerical  Results  6o 

Marshall-Palmer  Distribution  — 6l 

Polyakova-Shifrin  Distribution 63 

5.0  PERFORMANCE  OF  CANCELLATION  SYSTEM  WITH  UNIFORM  RAIN 

CHARACTERISTICS  68 

5.1  Rain  Characteristics — 68 

5-2  System  Equations  68 

5.3  Evaluation  of  Performance  — — 76 

5.1*  Comparison  of  Performance 80 

6.0  PERFORMANCE  OF  CANCELLATION  SYSTEM  WITH  A REALISTIC  KAIN 

MODEL 88 

<5.1  The  Rain  Model  — 88 

6.2  System  Equations 89 

6.3  Lower  Bound  on  Performance  — 99 

6.1*  Upper  Bound  on  Performance  — 108 


7.0  CONCLUSIONS 


REFERENCES 


— 113 


Appendix  A:  ANALYSIS  OF  MOMENTS  FOR  RANDOM  SUMS > — ll6 

Appendix  B:  THE  MEAN  VALUE  OF  Nk 121 

Appendix  C:  EVALUATION  OF  THE  FOURTH  MOMENT  OF  THE  RADIATION 

CHARACTERISTIC 123 

Appendix  D:  STATISTICAL  CONSIDERATION  OF  A PRODUCT  126 

Appendix  E:  THE  CROSS  CORRELATION  OF  X_(k),  X.(k),  Y.(k),  AND 

Y (k) ± i 128 

2 

Appendix  F:  THE  REDUCTION  IN  DIMENSION  OF  INTEGRALS  USED  TO 

EVALUATE  THE  EXPECTATION  OF  G'(k)  AND  |G'(k)|2  13U 

Appendix  G:  MOMENTS  OF  THE  S MATRIX  ELEMENTS 138 

Appendix  H:  NUMERICAL  INTEGRATION  TECHNIQUES 1^5 


».Mu4*,*>t  • *n?*'  ^ <*  *V»Mt*-4M *y^*  ** '.'  ? 


EVALUATOR'S  MEMO 


Project  No:  AS060188 


Contract  Mo:  F30602‘7k<-0037 


Effort  Title; 


Contractor : 


: Investigation  of  Radar 

Rain  Clutter  Using  a 
Polarization  Method 

University  of  Tennessee 
Knoxville,  Tennessee 


Radar  backseat ter  from  rain  rapidly  degrades  the  capability  to  detect 
aircraft  targets.  A most  crucial  situation  exists  with  GCA  (Ground 
Control  Approach)  and  PAR  (Precision  Approach  Radar)  equipment  during 
rain  storms.  This  effort  involved  a method  for  cancelling  the  backscatter 
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orientation  of  the  backscattered  polarization  and  in  turn  adjusts  the 
receiver  to  null  out  this  polarization.  The  target  return  will,  on  the 
average,  lose  about  3 d8  while  the  loss  in  rain  return  will  exceed  that 
from  ordinary  CP  by  about  10  dB. 

The  results  are  in  the  form  of  simulations  modeled  from  meteorological 
data.  With  a relatively  simple  scheme,  equipments  now  employing  CP 
can  be  converted  to  cancel  the  estimated  mean  polarization  of  the  rain. 
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1.0  INTRODUCTION 


1.1  Historical  Development 

In  spite  of  improved  performance  due  to  advances  in  detection 
methods,  low-noise  receivers,  waveform  optimization,  etc.,  the  modern 
radar  system  often  gives  a much  lower  level  of  performance  because 
of  rain  clutter. 

Early  radar  systems  vere  generally  designed  to  emit  and  receive 
linearly  polarized  radiation.  There  is  nothing  inherent  in  the  struc- 
ture of  the  linearly  polarized  system  to  suppress  the  clutter  caused 
by  precipitation.  Recognition  of  this  problem  prompted  investigation 
into  different  types  of  systems.  In  19^7  Ridenour  [lj  proposed  a 
significant  new  concept  for  the  basic  radar  system.  Ridenour  noted 
that  if  rain  particles  could  be  correctly  modeled  as  spheres  then  the 
use  of  circular  polarization  on  both  the  transmitting  and  receiving 
antennas  offered  the  possibility  of  significant  clutter  reduction. 

In  the  same  discussion  Ridenour  was  also  able  to  offer  experimental 
data  that  Indicated  that  clutter  reduction  on  the  order  of  26dB.  was 
possible  using  circular  polarization.  Additionally  Ridenour  was  able 
to  point  out  that  two  limitations  of  this  approach  to  clutter  cancel- 
lation are  rain  particles  which  are  not  true  spheres  and  radiation 
which  is  not  exactly  circular.  In  spite  of  the  advantages  offered 
by  the  use  of  circular  polarization  it  was  not  immediately  adopted. 

In  fact,  circular  polarization  was  rediscovered  in  195^  by  White  [2] 
while  at  the  same  time  Hunter  [3]  reemphasized  the  limitations  of  the 
technique  due  to  radiation  which  was  not  perfectly  circular. 


A necessary  step  in  the  process  of  correcting  the  errors  intro- 
duced by  non-circular  polarization  and  non-spherical  scatterers  was 
that  of  obtaining  an  accurate  model  of  the  scattering  particles  and  then 
obtaining  a solution  to  the  scattering  of  electromagnetic  waves  by 
an  object  whose  shape  is  that  of  the  model.  From  the  literature  one 
recognizes  that  interest  in  modeling  rain  particles  preceeded  the 
development  of  radar  technology  by  many  years.  A recent  paper  by 
Pruppacher  and  Beard  [3+1  offers  a review  of  the  efforts  along  this  line 
while  another  paper  by  Pruppacher  and  Fitter  [5]  claims  to  offer  the 

most  accurate  model  yet  available  for  rain.  In  simp : terms,  Pruppacher 

and  Pitter  claim  that  the  oblate  spheroid  is  a good  approximation  to 
the  true  shape  of  rain  particles  in  many  cases.  Additional  detail 
concerning  particle  shape  will  be  discussed  in  Chapter  3. 

The  scattering  of  electromagnetic  waves  by  bodies  of  various  shapes 
is  a problem  that  has  also  received  considerable  attention  for  many 
years.  Logan  [6]  notes  that  some  of  the  earliest  work  done  in  this 
area  involved  scattering  by  a sphere  and  was  done  by  Clebsch  in  l86l. 

In  the  same  article,  Logan  offers  an  excellent  outline  of  the  histor- 
ical development  of  the  theory  of  scattering  by  spheres.  Another 
paper  deserving  mention  at  this  point  is  a survey  of  scattering  liter- 
ature by  Corriher  and  Pyron  [7].  It  should  be  noted  however  that 
none  of  the  sources  referenced  by  the  two  previously  mentioned  works 
considers  the  problem  of  scattering  by  an  oblate  spheroid.  This 
problem  was  apparently  first  solved  by  Labrum  [8]  in  1952,  and  later 
Atlas,  Kerker,  and  nx+schfield  ;9l  offered  concurring  results. 
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1.2  Description  of  the  Proposed  Method 


This  report  is  an  analysis  of  the  system  depicted  in  figure  1.2-1. 
In  the  block  diagram,  the  transmitter  is  represented  in  the  left  hand 
portion  of  the  figure,  the  propagating  medium  and  scattering  elements 
are  shown  in  the  center  of  the  figure,  and  the  receiver  is  depicted 
in  the  right  hand  portion  of  the  diagram. 

The  features  which  distinguish  this  system  from  the  ordinary 
circular  system  are  indicated  in  the  block  diagram  as  elements  T,  R, 
and  G.  Elements  T and  R account  for  the  magnitude  and  phase  imper- 
fections which  prevent  the  ordinary  system  from  being  truly  circularly 
polarized.  G is  a dynamic  gain  element  which  has  been  added  in  the 
effor+  to  reduce  the  performance  limitations  of  the  ordinary  system 
caused  by  non-spherical  rain  and  non-circular  polarization.  A detailed 
description  of  G will  be  given  in  section  5. 


ca: 
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2.0  BACKGROUND  THEORY  AND  ASSUMPTIONS 


Prior  to  detailed  development  of  the  concept  which  gives  rain 
clutter  cancellation  it  is  helpful,  indeed  necessary,  to  first  set 
down  certain  definitions,  assumptions  and  preliminary  theory  appli- 
cable to  the  whole  work.  For  example,  to  avoid  confusion  we  shall 
use  the  word  "signal"  to  imply  either  a waveform  resulting  from  back- 
scatter  from  a target  or  waveforms  preceeding  target  or  rain  illumi- 
nation while  the  word  "clutter"  will  apply  to  waveforms  resulting 
from  rain  backscatter.  Other  definitions  are  made  clear  in  the  course 
of  development-.. 

We  make  the  following  basic  assumptions  throughout  the  effort. 

1.  Only  rain  clutter  is  present.  Thus,  noise,  ground  clutter, 
multipath  and  other  interference  sources  are  neglected. 

This  assumption  is  approximately  satisfied  for  a low-noise, 
narrow-beam  system  directed  several  beamwidths  above  the 
horizon. 

2.  Only  a single  target  exists  within  the  radar  beam. 

3.  The  target  behaves  as  a point  reflector,  that  is,  the  target 
size  is  small  relative  to  the  pulse  length  in  space. 

U.  There  are  no  multiple  reflections  between  rain  particles 
or  between  the  rain  and  target. 

5.  The  radar  is  a monostatic  system.  Generalization  to  the 
bistatic  system  1b  relatively  straightforward  and  some  of 
the  developments  are  formulated  in  such  a manner  prior  to 
specialization  to  the  monostatic  case. 


w^"*1  i^r*;p7w>^'  wpTTm'.wwwii 


6.  The  system  is  capable  of  transmitting  and  receiving  two 
linearly  polarized  waves  that  are  approximately  space- 
orthogonal. 

7.  The  antenna  patterns  have  narrow  main  lobes,  that  is,  beam- 
widths  are  on  the  order  of  a few  degrees  or  less  in  any  plane 
containing  the  boresight  axis. 

8.  Antenna  pattern  sidelobes  fall  off  with  angle  such  that  the 
dominant  contribution  to  the  rain  clutter  power  comes  from 
angles  sufficiently  near  the  direction  of  the  mainlobe  maxi- 
mum that  small-angle  approximations  may  be  made. 

2.2  Coordinate  Definitions 


Figure  2.2-1  defines  the  coordinate  system  which  interfaces  the 
radar  with  its  environment.  The  radar  is  at  the  origin  of  the  X,  Y,  Z 
system  and  radiates  in  the  Z (boreright)  direction.  Rain  particles 

are  assumed  to  fall,  on  the  average,  along  a line  parallel  to  0-x. 

/ 

This  line  is  defined  by  angles  C,  called  the  rain  heading,  and  ij>, 
called  the  rain  fall  angle.  With  these  definitions  the  average  wind 
vector  lies  in  the  0-x-y  plane;  it  also  lies  in  the  X-Z  plane  (as 
shown)  if  the  boresight  axi3  is  horizontal  and  the  wind  is  entirely 
horizontal.  In  general,  even  for  horizontal  wind,  the  wind  vector 
is  tilted  in  the  O-x-y  plane  when  the  boresight  axis  is  elevated. 

The  above  definitions  apply  to  the  fall  geometry  for  raindrops 
which  is  determined  by  the  wind.  Locations  of  the  drops  are  defined 
by  spherical  coordinates,  radar  centered.  The  coordinates  for  the 
ith  drop  are  , 8^  analogous  to  C»  ♦ • 
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2.3  Scattering  Theory 

It  is  desirable  to  be  able  to  write  an  equation  relating  the 
electric  field  incident  upon  and  reflected  by  a scattering  body  such  as 
a rain  particle.  The  first  step  in  the  process  of  obtaining  the  de- 
sired equation  is  the  specification  of  the  nature  of  the  scattering 
particle.  For  the  purposes  of  the  following  calculations  it  shall  be 
assumed  that  the  scattering  particle  is  a homogeneous  oblate  spheroid. 
Justification  for  the  assumed  shape  is  provided  in  section  3. 

N.R.  Labrum  [8]  has  provided  a fundamental  solution  to  the  problem 
of  electromagnetic  scattering  by  an  oblate  spheroid.  His  approach 
to  the  problem  is  based  upon  work  done  by  Stratton  [10],  which  has 
shown  that  where  the  ratio 

particle  radius 
radiation  wavelength 

is  small,  the  dipole  moment  induced  in  a sphere  is  the  same  for  both 
static  and  time  varying  incident  fields.  Labrum' s work  relates  the 
dipole  moment  induced  in  the  scattering  particle  to  the  incident  elec- 
tric field.  In  turn, the  reflected  electric  field  may  be  related  to 
the  induced  dipole  moment,  but  this  is  unnecessary  for  the  current 
work  since  the  final  equation  will  involve  only  ratios  of  this  param- 
eter which  associates  a reflected  field  with  an  incident  field.  Thus, 
the  equations  relating  the  reflected  electric  field  to  the  incident 
electric  field  will  be  of  the  form 
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E“  » SE+  (2.3-1) 

r - reflected  electric  field 
E+  ■ incident  electric  field 

but  S will  not  be  unitless.  As  in  Labrum’s  work,  S will  have  the 
units 


dipole  moment 
electric  field  intensity 


Farad  - meter 


2 


(2.3-2) 


Rice  and  Peebles  (11}  have  extended  Labruo'B  work  to  the  case 
of  circular  polarization  and  homogeneous  oblate  spheroids.  From  the 
standpoint  of  notation,  this  extension  is  simplified  if  the  field 
intensities  are  resolved  into  orthogonal  components  which  may  be  entered 
as  elements  in  a column  matrix.  Thus  the  following  shall  be  used. 


(2.3-3) 


For  the  particles  assumed.  Rice  and  Peebles  define  the  terms  in 
the  following  manner. 


su  * gJ£  cos2(0)  + gy  sin2(e) 
S12  * (g*  - gy)  cos(e)  eo»(t) 
321  " ~ ®y^  cos(e)  cos(*) 
S22  * ®x  cob2^^  + Ey  »in2(*) 


(2.3-*0 
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8WCJ  p(l-p2)3/2  (Cg-^Jb3 

**  ^ [Z&p*  - qv  + 2p arc tan  (p/A-p*))  Ug-e^  * 2e1(l-p2)3^2 


8**2  p(l-p2)3/2  (tg-e^b3 

^ ^ ^ - p/l-pJ  - arctan  (p//l-p2)]  (e~-«^7"+  2e1(l-p2)3^2 


p * a/b 

In  equation  (2.3-U),  0 1b  a particle  orientation  parameter  which 
is  used  by  Labrum  and  is  related  to  the  coordinate  parameters  used 
in  this  work  as  indicated  in  equation  (2.3-5). 


cos(6)  ■ sin(c)  sin(f) 


(2.3-5) 


The  results  of  applying  equation  (2.3-5)  to  the  definition  of  the 
S . . terms  given  in  equation  (2.3-10  are  provided  in  equation  (2.3-6) 


sn  * «y  ♦ («*  - Sy ) sin2  (Osin2  (♦) 


S12  * ^8x  “ 8,r)8in(c)ain(<»)cos(<>) 


(2.3-6) 


S22  " *x  co>2  ^ + *yBin2  W 


With  the  8 matrix  elements  defined  as  in  equation  (2.3-6),  it 
is  now  possible  to  express  the  components  of  the  scattered  electric 
field  in  terms  of  the  components  of  the  incident  electric  field  as 
indicated  in  equation  (2.3-3).  Note  that  the  equation  in  this  form 


J 


represents  scattering  from  a single  scattering  particle.  The  effects 
of  multiparticle  scattering  will  be  represented  by  a summation  of 
the  following  form. 


Nk  \ 

K - * (&!„  E . + ) + I (S,^  E /) 


i«l 


111  xi  ' ^ 12i  yi 


Nk  \ 

V * £ <S21i  O * £ (S22i  V+) 


(2.3-7) 
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In  equation  (2.3-7),  the  total  rnsnber  of  scattering  particles  within 

the  scattering  volume  of  interest  is  N . 

k 


2.k  The  Reference  System 


What  will  be  called  the  reference  system  in  this  work  is  presented 
in  block  diagram  form  in  Figure  2.1<-1.  The  system  iE  a linearly  polar- 
ized monostatic  radar,  and  the  performance  of  the  ordinary  circularly 
polarized  system  and  the  cancellation  system  will  be  compared  to  the 
performance  of  this  reference  system. 

From  the  theory  of  linear  antennas  [12],  it  can  be  shown  that  the 
electric  field  intensity  in  the  far  zone  of  a dipole  radiator,  e+, 
is  related  to  the  driving  point  potential,  e,p,  as  shown  in  equation 
(2.U-1). 


+ 


e 


, *1  eT  GT 
Z 


-J80Z 


(2.U-1) 


In  the  preceeding  equation,  is  a complex  constant,  6^  is  the  tra- 
ditional phase  constant,  Z 1b  the  radial  distance  between  the  phase 
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center  of  the  antenna  and  the  point  where  the  field  intensity  is  to 
he  measured,  and  is  the  radiation  characteristic  (voltage)  of  the 
transmitter  antenna. 

By  the  scattering  theory  of  the  preceeding,  section,  the  scattered 
field  at  a distance  Z'  from  the  scattering  body  has  an  intensity, 
e , which  is  given  by 


Sn e -jv 

e * e 


(2.4-2) 


where  is  an  element  of  the  previously  defined  scattering  matrix. 

Prom  the  theory  of  reciprocity  as  it  applies  to  linear  antennas, 

the  received  signal,  e , is 

r 


- K2GrSlle  -1* 0 2' 

e - K_G  e - e 0 

z- 

K2KlCrGTSlleT  -J60(Z+5r) 

* — — — - ■■  '■  e 

Z*Z" 


(2.4-3) 


For  a monostatic  radar,  Gy  ■ GT  * G and  Z » Z' , and  to  simplify  notation 
let 


K “ K2  ' K1 


(2.4-4) 


Thus  the  signal  received  as  a result  of  scattering  from  one  scat- 
tering body  is  given  by  equation  (2.4-5). 


~7~ 


Extending  the  analysis  to  cover  mult ipar tide  scattering  involves  only 
implementing  a summation  of  returns  from  the  individual  scattering 
bodies.  Such  a summation  is 


> Z 
i-1 


(k)gi 

r,2 


-J2B0Zi 


(2.U-6) 


The  receiver  signal,  e , indicated  in  equation  (2 .4-6)  is  a re- 

r 

suit  of  returns  from  all  the  scattering  bodies  lying  within  the  range 
cell  corresponding  to  the  time  at  which  er  is  evaluated.  A convenient 
way  of  noting  this  fact  is  to  specify  the  location  of  each  particle 
with  respect  to  the  center  of  the  range  cell  in  which  it  lies.  For 
"til 

the  i — particle  in  the  k — range  cell  this  would  be 


Zt(k)  » Z(k)  ♦ AZa  (2.14-7) 

An  approximation  which  will  be  used  in  the  denominator  of  equation 
(2.U-6)  is 

Z1(k)  s Z(k).  (2.U-8) 


This  approximation  is  in  error  by  no  more  than  0.5  percent  for  a trans- 
mitter pulsevidth  of  one  microsecond  and  a range  of  16  km.  Application 
of  the  approximation  to  equation  (2.U-6)  results  in  the  following. 

Note  that  is  the  total  number  of  scattering  particles  in  the  k^ 

range  cell. 
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(2.14-9) 
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In  sunnary,  equation  (2.4-9)  is  an  expression  for  the  output  voltage 

i.u 

resulting  from  clutter  oceuring  in  the  Yc~-  range  cell  for  what  will  be 
called  the  reference  system. 

2.5  The  Ordinary  System 


An  ordinary  system  will  be  defined  as  a radar  which  radiates  and 
receives  the  same-sense  circular  polarization  to  reduce  rain  clutter. 

Such  a system  is  ideal  if  it  radiates  perfectly  circular  polarization 
of  one  sense  and  has  zero  response  to  a perfectly  circular  wave  of  the 
opposite  sense.  An  imperfect  ordinary  system  does  not  produce  perfectly 
circular  operation  due  to  manufacturing  and/or  alignment  errors. 

By  recognizing  that  a circularly  polarized  wave,  indeed,  even  an 
elliptically  polarized  wave,  may  be  considered  to  be  comprised  of  two 
linearly  polarized  waves  we  may  model  the  transmission  and  reception 
operations  of  the  system  as  shown  in  Figure  2.5-1.  A signal  eT  is  power- 
split  with  one  signal  driving  an  "antenna"  representing  the  x-polarized 
pattern.  The  second  signed  passes  through  a complex  gain  constant  T, 
representing  the  error  that  the  transmitting  function  introduces  relative 
to  true  circular  polarization  (T*lt  ideally).  This  signal  then  passes 
through  a * ir/2  or  - tr/2  phase  shifts  to  produce  either  left  or  right 
hand  polarizations  respectively  when  transmitted  through  an  "antenna" 
representing  the  Y-polarized  pattern.  Similar,  but  reverse,  operations 
exist  in  the  receiver.  The  parameter  R of  the  receiver  equals  T in 
a monostatic  system. 

If  the  true  radiations  e*,  and  e*,  are  non-orthogonal , the  compo- 
nents each  couple  to  the  desired  space  axes  X and  Y (see  Figure  2.2-2). 

The  coupling  may  be  accounted  for  by  use  of  a coupling  matrix  with  elements 


t ^ , i and  and  2.  Thus, 
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(2.5-1) 


Using  Figure  2.2-2  we  may  find  that 
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*11 

*12 

cos(ex) 

-sin(6y) 

*21 
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sin(0x) 

cos(0  ) 

y J 

(2.5-2) 


Another  matrix  accounts  for  either  rain  or  target  scattering. 

Still  a third  matrix  accounts  for  the  receiver  non-orthogonal  polariza- 
tions. For  a monostatic  radar  the  matrix,  defined  to  have  elements 
r^j,  i and  J=1  and  2,  relates  received  fields  at  the  antenna  to  reflec- 
tions from  range  cell  k as  follows : 


v(k) 

ii 

V11 

r!2 

e”(k) 

e Y- (k) 

_r21 

r22 

e”(k) 

where , 

cos(°x) 

-sin(0y) 

cos(ex-0Y^ 

Here  [•]*  represents  matrix  transpose. 

We  now  calculate  the  received  waveforms  for  the  system  of  Figure 
2.5-1.  We  shall  use  the  exponential  (phasor)  representation  of  wave- 
forms. The  transmitter  excitation  is  then 
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where  1^  is  the  magnitude  of  e^  and  uig  is  radian  frequency*  Let 
G+(a,  3)  and  G^o,  3)  he  the  transmit  patterns  for  X'  and  Y"  polari- 
zation components  evaluated  in  directions  o and  3.  The  electric  fields 
e+"  and  e^  near  the  antenna  become 

V " 0)  (2.5-6a) 

V « ♦JTK^a,  3)  eT/^'  (2.5-6b) 

where  1^  is  a proportional ity  constant. 

To  continue  further  we  first  assume  scattering  is  from  a single 
rain  particle.  Results  are  then  extended  to  multiple-particle  rain 
and  target  scattering. 

Let  rain  particle  i be  located  in  range  cell  k at  a distance 


Zt(k)  - Z(k)  + AZt 


(2.5-7) 


where  Z(k)  is  the  distance  to  the  center  of  the  cell  and  AZ^  is  particle 
position  relative  to  cell  center.  By  combining  (2.5-5)  with  (2.5-1) 
the  field  components  e*  and  e*  may  be  determined.  If  we  then  use  the 
scattering  matrix  of  the  rain  particle  and  include  the  attenuation  fac- 
tor due  to  dist<  - ce  we  may  write  the  received  fields  e”(k)  and  e.^(k) 
due  to  range  cell  k as 
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(2.5-6) 
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Here  the  rain  scattering  matrix  [a]  has  elements  sinn£ (lc)  for  cell  k. 

We  may  next  use  (2.5-3)  with  2.5-8)  to  obtain  the  fields  in  the 
polarization  planes  of  v,he  receiving  antenna.  The  received  waveforms 
e^(k)  and  e2(k)  are  related  to  these  fields,  e^-(k)  and  e“4k),  by 


e1(k)  **  KgG*(«^ , 3^  e^'(k) 


(2.5-9a) 


(2.5-9b) 


e2(k)  = +JKgRG*(ai,  Bi)  e">(k) 

where  Kg  is  a proportionality  constant . Combining  and  recognizing 
XXX  Y Y Y 

that  R=T,  G . = G = G and  G . = G = G for  a monostatic  system  gives 
* t r t r 
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(2.5-10) 


where  K * 1^  • Kg. 

To  be  more  explicit,  (2.5-10)  may  be  expanded  to  obtain  e^(k)  and 


e?(k)  as 
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Certain  special  cases  of  these  results  are  of  interest.  For  example, 
most  of  our  work  will  Involve  assuming  the  system  polarizations  are 
orthogonal  (8^  = 0 and  6^  «*  0 in  Figure  2.2-2);  with  this  assumption 

*11  * *22  “ 1 and  *12  * *21  " °*  Hence, 


e1(k)  * GX(u1,  Pj)  tsnl(k)  GX(ai , 81>  +JTs12i(k)  GY(ai,  B± ) } 
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Another  special  case  derives  from  an  assumption  that  the  patterns 
X y 

G and  G are  identical.  We  have 


GX(a1,  Bj)  - GY(a1,  Bj)  - G^,  6^ 


(2.5-13) 


giving 
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(2.5-14b) 


Any  of  the  single  particle  rain  clutter  results  may  be  used  for 
the  multiple  particle  case  by  simply  summing  over  all  particles  in 
cell  k.  The  result  of  such  a summation  coupled  with  the  simplification 
indicated  in  equation  (2.U-8)  is  given  below  for  equation  (2.5-12). 
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Similarly,  (2.5-1**)  becomes 
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In  all  these  target  signals  ve  shall  most  often  assume  the  target  is 
located  on  the  beam  axis  so  X « 0 and  Y ■ 0. 


1 4 
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3.0  THEORETICAL  RAIS  MODEL 


In  this  section  a rain  model  is  developed  which  is  used  in  sec- 
tions U,  5 and  6 for  calculating  performance  of  ordinary  and  cancellation 
systems.  The  model  is  built  around  the  specification  of  raindrop  shape, 
size  distribution  and  orientation  with  respect  to  wind  velocity  com- 
ponents . 


3.1  Drop  Shape 

It  is  well  known  that  falling  raindrops  are  not  spherical  [U,  l'U, 

16,  20 ].  The  drops  become  flattened  on  the  bottom  as  they  fall  and  may 
even  develop  a small  "dimple"  at  the  center  of  the  flattened  surface  [13] 
as  shown  in  Figure  3.1-l(a).  Very  small  droplets  fall  with  the  least 
distortion.  Distortion  becomes  more  severe  the  larger  the  raindrop 
and  is  significant  for  drops  of  2uim  equivalent  spherical  diameter^- 
or  larger.  Above  about  2 mm  the  dimple  occurs  and  becomes  more  severe 
as  equivalent  spherical  diameter  increases  up  to  about  U.5mm  where 
drops  tend  to  break  up,  even  in  quiet  air,  [15). 

Being  more  specific,  it  is  found  [5]  that  for  equivalent  spherical 
diameters  less  than  about  0.17mm  the  drop  is  spherical  with  negligible 
error.  Between  0.17nm  and  0.5nm  the  drops  are  oblate  spheroids.  Above 
Q.5mn  to  about  U.5  cm  the  drops  are  asymetrical  spheroids. 


1 


The  equivalent  sphere  is  that  having  the  same  volume  as  the  distorted 
drop. 
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Pruppacher  and  Pitter  [5]  have  given  an  excellent  model  for  rain- 
drop shape.  It  appears  to  fit  measured  results  quite  veil  except  for 
the  largest  raindrops  vhere  the  actual  deformations  are  slightly  under- 
estimated hy  the  model.  We  shall  next  describe  the  shape  model  of 
Pruppacher  and  Pitter  (P&P  model)  and  Justify  that  our  use  of  the  ob- 
late spheroid  is  a reasonable  second-order  approximation  to  the  P&P 
model. 

The  P&P  model  is  a refinement  of  work  due  to  Savic  [17].  Actual 
drops  will  have  a cross-section  such  as  shown  in  Figure  3.1-l(a). 

The  drop,  of  coursek  has  rotational  symmetry  about  the  vertical  axis. 

The  radial  distance  to  the  surface  from  some  interior  point  is  r(0). 

The  P&P  model  consists  of  expanding  r(8)  as  an  arbitrary  (but  neces- 
sarily single- valued)  function  into  a Fourier  series1 


r (0 ) " b0[l  + ” 
n«0 


c 

n 


cos(n0)] 


(3.1-1) 


and  finding  the  "distortion  coefficients"  cn  which  correspond  to  a 

balance  of  forces  (aerodynamic  pressure,  hydrostatic  pressure  gradients 

In  the  drop,  pressure  across  the  surface  and  surface  tension)  on  the 

falling  drop.  Here  bQ  is  the  radius  of  the  equivalent  sphere  having 

the  same  volume  as  the  deformed  drop  and  the  c were  found  under  a 

n 

constant-volume  constraint.  It  results  that  c^  * 0 may  be  assumed 
without  loss  in  generality. 
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Notation  of  [5]  has  been  altered  to  agree  with  that  used  here. 
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The  solutions  of  [5]  for  the  cr  values  are  plotted  in  Figure  3*1-2 
as  a function  of  b^.  The  actual  deformation  ratio  b/a,  with  b and  a 
defined  in  Figure  3.1-l(a),  is  a function  of  b^  as  shovn  in  Figure  3.1-3* 
Notice  how  excellent  the  fit  is  to  experimental  data. 

We  now  show  that  an  oblate  spheroid  is  a reasonable  approximation 
to  the  more  exact  P4P  model.  The  reason  for  using  the  oblate  spheroid 
in  analysis  is  that  analytical  solutions  for  the  backscattered  fields 
exist  [8]  in  relatively  simple  forms. 

Consider  expanding  r(6)  for  an  oblate  spheroid  having  the  cross- 
section  of  Figure  3.1-l(b).  Since  the  cross-section  is  an  ellipse 


r2(e) 
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e 


2 

sin2(0) 


(3.1-2) 


where 


e2  - 1 - (b/a)2  < 1 


(3.1-3) 


is  the  eccentricity  of  the  ellipse.  We  select  the  two  free  parameters, 
a and  bt  for  the  oblate  spheroid  first  to  make  its  volume  equal  to  that 
in  the  P&P  model  and  second  to  make  the  ratio  b/a  nearly  equal  to  that 
of  the  P4P  model.  In  the  latter  case,  the  actual  relationship  between 
b/a  and  b^  is  the  solid  line  cf  Figure  3*1-3.  We  shell  choose  b/a 
according  to 


which  is  shown  dashed.  This  function  is  a good  approximation  to  the 
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Fig.  3.1-3.  Variation  of  drop  deformation  with  drop  size: 
Comparison  of  theory  and  experiment.  Figure  from  Pruppacher 
and  Pltter  [5j. 
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P&P  model  except  for  large  where  it  actually  approximates  the  measured 
data  better  than  the  P&P  model. 

By  equating  the  volume  kira^b/3  [ 18 ] of  an  oblate  spheroid  to  the 
volume  hirbQ  /3  of  an  equivalent  sphere  we  find  that 

b3  - b3(l  - e2)  (3.1-5) 


for  a constant  volume  constraint  and  any  e < 1.  The  constraint  that 
e must  be  chosen  to  approximate  that  of  the  P&P  model  requires  that 
e be  given  by  (3.1-3)  using  (3.1-*0.  From  (3.1-2)  using  (3.1-5)  we  get 


a »g  (1  - s2>2/3 

r <e>  

1 - e^sin  (0) 


(3.1-6) 


This  function  has  even  symnetry  and  has  two  cycles  of  behavior  in  the 
interval  (-it,  it).  Thus,  the  Fourier  series  representing  (3.1-6)  has 
only  terms  for  n * 0,  2,  U,  . . . when  the  expansion  is  over  (-tt,  it). 
The  series  may  be  written  in  the  fern 


r(0)  * b.  (l  + cA  + E c cos(n0)} 
0 0 _ n 

n*2 
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(3.1-7) 


Solutions  for  the  coefficients  cn  involve  elliptic  integrals  and 
are  slightly  involved  to  determine.  The  first  two  are 
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where  K(e)  and  D(c)  are  complete  elliptic  integrals  defined  by  C 19 1 
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(3.1-9) 


(3.1-10) 


We  may  now  demonstrate  that  the  oblate  spheroid  is  a reasonable 
approximation  to  the  PtP  drop  shape  model.  By  comparing  (3.1-7)  for 
the  oblate  spheroid  with  (3.1-1)  for  theP&P  model  we  see  that  the  cr 
should  be  nearly  equal  to  the  c^  for  all  n.  Two  values  of  c^,  as  given 
in  (3.1-8),  are  plotted  in  Figure  3.1-2  as  dashed  curves.  Excellent 
agreement  is  seen  for  all  even  values  of  n.  The  main  loss  in  using  the 
oblate  spheroid  is  the  loss  of  odd  numbered  terms  for  n >_  3.  There 
is  only  one  odd  term  in  the  PfcP  model  of  any  significance  (n  * 3)  how- 
ever, so  the  loss  is  small.  (The  term  for  N * 5 is  well  over  two  orders 
of  magnitude  smaller  than  unity  for  all  values  of  b^). 


3.2  Drop-Size  Distributions 


The  various  mechanisms  that  cause  rain  to  have  a given  drop-size 
distribution  are  so  complicated  as  to  prohibit  analysis.  Indeed,  there 
is  such  variation  in  measured  distributions  in  the  literature  that 
one  may  conclude  that  no  one  distribution  can  describe  all  rain  cases. 

The  best  that  one  can  hope  for  is  to  find  a raindrop  distribution  model 
which  is  simple,  easy  to  work  with,  and  represents  a more  or  less  average 
of  many  rain  situations.  With  these  facts  in  mind  we  shall  discuss 
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various  measured  data  from  the  literature,  consider  available  proposed 
models,  and  finally  select  the  model  which  seems  to  present  the  best 
compromise  in  the  desired  behavior. 

Before  proceeding  it  is  helpful  to  classify  the  forms  of  distri- 
butions according  the  Figure  3.2-1.  We  define  these  forms  as  types  A, 

B,  C,  and  D.  The  general  behevior  is  all  that  is  important  in  these 
definitions  and  some  variations  may  be  expected. 

Measured  Data 

Measurements  of  raindrop  sizes  were  first  made  as  early  as  1895 
in  Germany  by  J.  Wiesner  using  an  absorbent  paper  method.  For  our 
purposes,  however,  the  most  significant  early  paper  was  due  to  Laws 
and  Parsons  [21]  who  used  a flour  pellet  method  developed  by  W.A.  Bentley 
in  190U,  an  American.  Their  data,  taken  in  1938  and  1939  in  Washington, 
D.C.  on  surface  rain,  were  approximately  type  C below  drop  diameters 
of  about  1 (at  2.5mm/hr  rain  rate)  and  2 (25mm/hr)  mm,  and  type  B above 
these  diameters. 

Marshall  and  Palmer  [22]  analyzed  surface  data  taken  in  Ottawa, 
Canada  [23]  in  19^6  using  the  filter  paper  method.  They  obtained  a good 
fit  to  their  data  and  data  of  [2l],  for  1 <_  I <_  25  mm/hr  using  their 
famous  distribute ci 

"Xd0  -3 

n(dg)  = Ae  m mm  (3.2-1) 

where 

A * 8000m”  3mm  ^ 


(3.2-2) 
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X - k.l  r0,21™"1  . (3.2-3) 

•J 

Here  n( ) ig  the  number  of  drops  per  unit  volume  (in  nr ) having  diameters 

from  dg  to  d^  + dd^,  I is  rain  intensity  in  mm/hr  of  water  rate  and 

dg  is  in  mm. 

Blanchard  [2k]  t using  the  filter  paper  method,  has  given  average 
distributions  for  a large  number  (113)  of  measurements  taken  in  Hawaii 
in  1952.  Some  of  the  data  were  taken  et  or  near  the  dissipating  edge 
of  non- freezing  orographic  clouds.  Other  data  corresponds  to  rain  at 
the  base  of  a non-freezing  orographic  cloud  while  some  applied  to  non- 
orographic  rain  (probably  due  to  snow  falling  through  the  melting  level). 
Orographic  data  were  mainly  type  B for  low  rain  rates  (<2mm/hr)  and 
type  D for  higher  rates.  The  non-orographic  data  were  type  C. 

Mason  and  Ramanadham  [25]  have  given  surface  data  on  rain  taken 

apparently  in  England  with  a photolelectric  spectrometer  which  is  quite 
similar  to  that  in  [UJ. 

Monsoon  rain  data  in  India  (1956)  have  been  reported  by  Murty  and 
Gupta  [26]  which  resulted  from  the  filter  paper  method.  Average  drop- 
size  distribution  curves  based  upon  229  sets  of  measurements  were  given 
for  rain  intensities  I from  2 to  1^0  mm/hr.  The  data  were  type  C and 
quite  similar  to  the  non-orographic  data  of  [2^]. 

One  of  the  most  comprehensive  studies  of  raindrop  distributions 
is  discussed  by  Fujiwara  [27,  28].  Data  taken  with  a raindrop  camera 
related  to  data  taken  at  Miami,  Florida  and  Champaigne,  Illinois  on 
a large  number  of  storms.  The  data  are  too  numerous  to  discuss  in  de- 
tail. However,  one  of  the  most  important  observations  to  be  made  is 
on  how  variable  drop-size  distribution  can  be.  It  can  vary  from  u 
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unimodal  structure  ( vlth  maximum  at  about  1 to  2 mm  diameter)  to  multi- 
modal with  many  peaks;  of  course,  the  exponential  decreasing  structure 
(type  B)  is  also  found.  The  variation  (fluctuations)  may  even  be  pro- 
nounced for  a single  distribution  measurement  so  results  do  not  usually 
form  smooth  functions.  Another  important  observation  is  that  distributions 
always  decreased  for  small  drop  sizes  (type  A).  Fuji war a also  shows 
time  evolution  of  distributions  which  indicate  great  variability  in 
only  short  time  periods  (as  small  as  1 min.).  Most  of  the  data  of 
[27]  were  for  thunderstorms.,  rainshowers  and  continuous  rains. 

More  data  on  tropical  rains  are  given  by  Sivaramakrishnan  [29]. 

The  filter  paper  mtehod  was  used  to  measure  surface  drop  distributions 
at  Poona,  India.  The  data,  for  non-freezing  rain  and  rain  from  melting 
snow  were  generally  of  type  B hut  considerable  fluctuation  existed  along 
the  curves. 

Additional  data  on  rain  from  a melting  band  is  given  by  Hardy  [30], 
Relatively  steady  rain  at  Flagstaff,  Arizona  on  31  July  1961  produced 
distributions  of  type  A for  I from  0.11  to  3.4  nun/hr.  These  results 
were  mostly  smooth  curves  with  the  peak  occuring  near  1 mm  drop  diameter. 
Measurements  were  made  with  a photoelectric  raindrop-size  spectrometer. 

The  photoelectric  spectrometer  was  also  used  by  Dingle  [31]  to 
measure  thunderstorm  rain  of  8 October  1959  at  Ann  Arbor , Michigan. 

Results  were  type  B where  the  slope  decreases  as  I increases  for  I from 
0.18  to  51. b mo/hr.  A Marshall-Palmer  (M-P)  model  is  a farily  good 

fit  to  these  data.  Other  data  for  a heavy  shower  on  23  October  1959 
and  for  showers  of  16  June  1960  also  indicated  reasonable  agreement  with 
the  M-P  distribution,  especially  when  the  data  were  averaged  (16  June 
results ) . 
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Using  at  3.2  am  doppler  radar  during  i960  and  1961  at  Pershore, 
England*  Caton  [32]  analyzed  83  of  107  observations  of  continuous  rain 
from  melting  levels  in  an  excellent  paper.  The  rains  were  related  to 
three  warm  fronts,  four  cold  fronts  and  three  low  pressure  centers. 
Melting  level  varied  from  1.0  to  3.3  km  above  ground.  Humidity  was 
nearly  always  high  and  evaporation  below  750  m was  negligible.  Median 
distributions  at  altitude  are  quite  smooth  in  their  behavior,  which  is 
type  A.  The  M-P  distribution  fits  the  data  only  over  the  mid-range  of 
drop  diameter  and  overestimates  both  small  and  large  diameter  drop  num- 
bers. Caton  also  gives  bounding  curves,  within  which,  905?  of  all  drop 
concentrations  fall.  These  bounds  are  also  smooth  curves.  Data  were 
for  I from  0.1  to  5.6  mm/hr. 

Distributions  for  high-intensity  rains  have  been  assembled  by 
Blanchard  and  Spencer  [33].  Their  data  are  taken  from  Mueller  and  Sims 
[3^]  and  Hudson  [35].  Data  of  [3*0  were  taken  from  several  types  of 
rain  using  a raindrop  camera  at  several  locations  around  the  world 
(Miami,  Florida;  Franklin,  North  Carolina;  and  the  Marshall  Islands). 

One  set  of  three  curves  representing  about  10  sample  distributions 
(type  A),  and  corresponding  to  I of  170,  205  and  216  ram/hr,  show  close 
agreement  to  the  M-P  function  for  dorp  diameters  from  1 to  about  3.5  nr.. 
For  diameters  outside  this  range  the  M-P  function  overestimates  the  drop 
density.  A curve,  for  I = 300  ram/hr,  of  data  from  [35]  were  quite 
similar  to  the  data  near  200  mm/hr  from  [?U]  (again  type  A,  the  curve 
averaged  9 samples).  For  I = 500  mra/hr,  data  based  on  averaging  3 
samples  [35]  again  gave  similar  performance  except  the  peak  moved  from 
near  1 mm  diameter  out  to  near  2 run.  All  data  from  [35]  were  taken 
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at  Miami,  Florida. 
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Discussion  of  Measurement  Data 

Including  Russian  data,  to  be  commented  upon  subsequently,  we  have 
discussed  measurements  from  all  >ver  the  world  (U.S.,  England,  Canada, 
Hawaii,  India,  Russia,  Africa,  *wad  the  Marshall  Islands)  covering  inten- 
sities from  a few  tenths  of  a am/hr  to  500  ana/hr.  The  data  apply  to  rrin 
at  various  spatial  positions  (within  rain  cloud,  at  rain  cloud  base, 
at  edge  of  cloud,  between  cloud  base  and  earth's  surface  and  at  the 
surface)  and  correspond  to  nearly  every  type  of  rain  generation  mech- 
anism ( from  melting  band  for  snow,  hail,  sleet,  fror  non-freezing  clouds, 
frontal  stormB,  showers,  convective  storms  and  both  freezing  and  non- 
freezing  orographic  rain).  Thus,  it  can  be  concluded  that  sufficient 
results  are  given  to  represent  nearly  all  typical  rain  situations  which 
can  occur. 

The  factors  which  cause  a particular  rain  situation  cure  extremely 
complex.  Seme  of  these  are:  temperature  and  humidity  profiles,  time 

and  spatial  location  of  points  of  Interest,  spatial  extent  of  rain  cloud, 
intensity,  accretion  (important  in  non-freezing  clouds),  coalescence  and 
drop  breakup  (related  to  intensity),  evaporation  (related  to  humidity), 
orography,  and  the  generation  mechanism.  Because  of  the  complex  nature 
of  the  problem  we  make  no  effort  to  determine  why,  where  or  how  a par- 
ticular situation  arises.  We  only  here  attempt  to  correlate  the  best 
rain  model  to  fit  an  average  rain  situation.  Thus,  we  place  the  averaged 
data  of  [20],  [24],  [26],  130],  [31],  [32]  and  [33]. 

It  is  helpful  to  firBt,  divide  all  caseB  into  two  categories : 
orographic  and  non-orograpbic.  Typical  distributions  in  the  former  case 
are  markedly  different  from  those  of  the  latter.  They  typically  behave 
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as  type  B (with  steep  slope)  for  low  Intensity  rain  and  as  type  P for 


heavier  rains. 


For  non-orographic  rain  it  is  again  convenient  to  place  data  in 


two  categories.  The  first  contains  restilts  obtained  by  the  flour  pellet 


and  filter  paper  methods.  The  second  has  data  from  photoelectric  spec- 


trometer, radar  and  raindrop  camera  methods.  An  immediately  obvious 


difference  in  distribution  form  between  the  two  categories  is  seen. 


For  small  drop  diameters  the  former  shows  increasing  concentrations 


(type  C)  while  the  latter  shows  decreasing  concentraxions  (type  A). 


One  is  tempted  to  suspect  that  splashing  of  large  drops  causes  the 


increase  in  the  filter  paper  method  (perhaps  even  the  flour  method). 


However,  as  pointed  out  by  Caton  [32]  the  splashing  problem  has  been 


considered  by  at  least  one  researcher  and  was  ruled  out.  It  is,  of 


course,  possible  that  both  categories  of  measurement  method  are  correct. 


On  the  other  hand,  a radar  man  tends  to  place  more  confidence  in  radar 


data.  Caton  has  carefully  considered  his  radar  system  and  concludes 


the  weight  of  evidence  supports  distributions  of  the  second  category 


for  moderate  and  heavy  rains.  Regardless  of  how  the  above  suspicion 


is  resolved,  at  least  one  thing  is  clear  from  nearly  all  of  the  data 


of  both  categories:  below  some  drop  size,  typically  in  the  1 to  3 mm 


diameter  range,  actual  rain  drop  densities  tend  to  be  overestimated 


by  the  M-P  distribution  (An  exception  would  be  very  light  rain  during 


very  high  humidity). 


In  the  mid-range  of  raindrop  size  the  M-P  distribution  is  a fair 


to  good  approximation  to  all  measurement  data.  The  region  of  close 


agreement  varies  with  intensity  and  Dyer  [36]  suggests  it  is  from  0.75 


to  2.25  mm  for  near  1 irzn/hr,  1.25  to  3 mm  for  near  5 mm/hr  and  1.5  to 


fc.5  non  for  near  25  mm/hr.  In  the  light  of  data  shown  heret  those  uppc  ' 
limits  are  perhaps  optimistic. 

For  larger  drops,  measured  distributions  almost  always  are  over- 
estimated by  the  M-P  curve.  This  fact  is  not  surprising  since  there 
is  a practical  size  limit  to  drops  caused  by  break-up. 

The  above  results  may  broadly  be  summarized  as  follows.  A reasonable 
model  for  the  drop-size  distribution  of  average  rainfall  has  these  char- 
acteristics: (l)  a peaked  behavior  occurs  with  the  peak  typically  oc- 

cur ing  to  0.5  to  1.0  nan  diameter  for  intensities  up  to  about  300  mm/hr, 
the  peak  moves  nearer  to  2.0  for  higher  intensities;  (2)  for  drop  di- 
ameters below  the  peak  region,  the  M-P  curve  overestimates  the  number 
of  drops;  (3)  for  drop  diameters  above  the  peak  region  and  up  to  a 
value  from  about  2 to  U nm  depending  on  intensity  (up  co  300  nm/hr) 
the  M-P  distribution  gives  results  nearly  equal  to  measured  valuer 
(typically  within  a factor  of  2);  (1*)  for  diameters  above  the  2 to  lumn 
region  the  actual  results  are  overestimated  by  the  M-P  curve. 

Prop-81 ze  Distribution  Models 

Having  determined  the  desirable  characteristics  of  a model  we  look 
at  the  various  models  suggested  in  the  literature. 

The  Marshall-Palmer  [22]  model  already  introduced  is  the  oldest 
and  still  remains  widely  used  because  of  its  simplicity  in  calculations. 
Another  reason  is  that  it  is  a good  approximation  in  the  mid-range  of 

drop  sizes,  which  contribute  greatly  to  radar  cross-section  because  of 
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the  dp  power  relationship.  An  often  played  game  Is  to  find  the  best 
A and  X to  fit  given  data.  Sene  examples  are  given  in  [36]  from  various 
researchers.  Perhaps  one  of  these  examples  should  be  mentioned.  Using 
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data  of  [26]  for  orographic  rain 


A = 7500  I 
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(3.2-U) 


A * 1*.5  1‘ 
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(3.2-5) 


are  derived.  Since  our  desired  model  characteristics  are  referred  to 
the  M-P  curve  ve  use  it  tc*  cor  pare  other  models. 

Best  [17]  gave  a rain  model  in  1950.  It  may  be  written  in  the  form 
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where 


a - AI* 


(3.2-7) 


W « Cl 


(3.2-8) 


with  A,  C,  p,  r and  n constants.  Here  W is  the  amount  of  liquid  water 
per  unit  volume  of  air  and  I is  rain  intensity.  If  <Iq  is  in  mm,  I in 
mm/hr  and  W in  mm  /nr  the  mean  constant  values  are  A * 1.3,  C = 67, 
p * 0.232,  r * 0.8U6  and  n * 2.25,  and  n(dQ)  is  in  m ran  . The  con- 
stants may  vary  appreciably  from  their  means,  especially  n if  rain  is 
showery  or  orographic  [37l. 

Taking  the  average  constant  values  (3.2-6)  becomes 
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Plots  of  (3.2-9)  show  clearly  that 
requirements  (it  is  more  like  type 
its  more  general  form,  however,  we 
simple  model  if  possible. 

In  1956  Litvinov  [38]  proposed 
in  [391*. 

“Ad0/2 

n(dQ)  ■ Me 

where  N and  X are  constants 
with  rain  type.  Thus,  for  a given 
since  it  does  not  depend  upon  I. 

Litvinov  [391  also  describes 
K.S.  Shifrin.  The  model  (P-S)  is 
[Ul].  It  is 

2 “Ydn 
n(d0)  - Ad2  e 


the  Best  model  does  not  fit  our  model 
C behavior).  It  may  be  better  in 
seek  not  only  a good  but  also  a 

the  following  model  as  discussed 


(3.2-10) 

as  far  as  intensity  is  concerned,  but  vary 

rain  type,  this  model  is  not  flexible 

«• 

It  is  not  considered  further, 
a model  due  to  E.A.  Polyakova  and 
also  described  by  Krasyuk,  et.  al. 

(3.2-11) 


where  A and  y are  given  in  Table  3.2-1.  Litvinov 
Table  3.2-1 


Type  Bain 

A (m~  8ram~  8 ) 

Y(n3n“1) 

Thawing  of  Pellets  (Hail) 

6U500I-°*5 

g.951-0.27 

Thawing  of  Granular  Snow  (Sleet) 

117501-°* 29 

U.87r0,2 

Thawing  of  MonGranular  Snow  (Snow) 

28201"°* 18 

i*.oir°*19 

This  article  was  translated  by  Mr.  Jim  Elliott  of  the  University 
of  Tennessee. 

**  la  a later  paper  by  Litvinov  [U0],  ar  translated  by  Irving  Bain, 

U and  X are  said  to  be  functions  of  intensity  but  not  of  the  usual 
form  al~®.  This  esse  would  appear  to  be  true  with  [30]  having  either 
a typographical  or  translation  error.  At  any  rate  (3.2-10)  does  not 
account  for  the  decrease  of  distribution  for  very  small  drop  size 
and  is  not  considered  further. 
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Fig.  3.2-5.  P-S  model  compared  to  measured  data  near  200  ram/h  rain 
intensity. 
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has  shown  Russian  data  for  all  three  types  of  rain  which  are  fit  excel- 
led ly  using  (3.2-11).  Figures  3.2-2,  3,  4 and  5 illustrate  the  P-S 
model  for  I * 2.5,  25,  tO  and  200  mm/hr  respectively.  Spotted  on  these 
figures  are  various  average  data  described  earlier;  also  shown  where 
possible  are  the  estimated  origins  of  the  rain,  that  is,  whether  orig- 
inating  from  melted  snow,  sleet  or  hail,  etc.  For  moderate  rain  , 
using  Figure  3.2-2,  there  is  excellent  agreement  with  averaged  data, 
especially  if  we  assume  the  Caton  data  (taken  from  August  to  December 
1961  in  England)  originated  as  hail.  For  heavy  rain  the  model  does 
not  fit  average  data  too  well  as  seen  from  Figure  3.2-3,  if  Dingle's 
data  are  from  melting  snow.  On  the  other  hand,  if  the  data  are  from 
melting  hail,  the  agreement  is  quite  good.  Similar  comments  may  be 
made  about  excessive  rain  from  Figure  3 2-!+.  For  deluge  rain  Figure 
3.2-5  shows  the  model  is  a poor  fit  regardless  of  the  rain's  origin. 

To  better  place  bounds  upon  the  usefulness  of  (3.2-11)  we  have  re- 
plotted some  of  Litvinov's  measured  data  [39]  in  Figures  3.2-6,  7 The 
f'irst  figure  is  for  rain  originating  from  hail.  Excellent  agreement 
is  seen  for  intensities  from  1.5  to  15.2  mm/hr.  (Litvinov  also  shows 
curves  for  3.47  and  7.25  ran/hr  which  also  show  close  agreement).  The 
curve  for  I - 1.5  represents  in  average  of  77  measurements  while  only 
a single  sample  applies  to  the  T ~ 15*2  curve.  The  second  figure  shows 
excellent  agreement  of  data  and  theory  for  rain  from  melting  snow  for 
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We  define  intensities  as:  drizzle,  0-0.5  mm/hr;  light  rain,  0.5— 2.0  ram/hr; 

moderate  rain,  2. 0-8.0  mm/hr;  heavy  r»+n,  8-32  nin/hr;  excessive  rain, 

32-128  mm/hr;  deluge  rain,  128  mra/h  higher. 
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Intensities  front  1.5  to  11.1  nan/hr  (curves  are  average  of  h and  ll» 
samples  respectively).  Excellent  agreement  is  obtained  for  intermediate 
intensities  of  3.52  and  7.^7  mm/hr  as  well  [39].  Finally,  Litvinov  also 
gives  results  for  rain  from  granular  snow  (sleet)  which  verify  agreement 
of  the  model  with  measurements  for  intensities  from  1.56  to  10.88  mm/hr. 

Thus,  in  summary,  the  P-S  model  appears  to  represent  average  rain- 
drop size  distributions  quite  well  up  to  intensities  of  at  least  15  mm/hr. 
There  is  some  doubt  as  to  its  validity  for  intensities  near  25  to  50  nm/hr 
end  it  is  not  a good  model  for  intensities  near  200  mm/hr.  It  can  prob- 
ably serve  useful  analytical  purpose  for  rains  through  the  heavy  region. 

Fujiwara  [27]  has  given  a model  involving  four  variable  parameters 
and  is  applicable  to  fitting  individual  storms.  It  is 
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(3.2-12) 


where  N is  the  total  number  of  drops  per  unit  volume;  DQ  is  a "mode 
shift"  parameter,  typically  around  0.8  mm  on  the  average;  S is  a "skew- 
ness" parameter  usually  varying  from  1 to  3;  a is  a "broadness"  param- 
eter which  is  related  to  rainfall  rate,  and,  from  fitting  of  measure- 
ment data,  ranges  from  1 to  3.5  nun.  Because  of  the  number  of  parameters 
having  unclear  relationships  with  intensity  and  rain  type  (3.2-12)  is 
difficult  to  work  with.  Fujiwara  shows  a good  fit  to  many  data  and 
(3.2-12)  may  be  very  useful  model  when  these  relationships  are  developed. 
However,  until  the  developments  are  available  we  shall  not  consider 
(3.2-12)  a useful  simple  model. 

Finally,  we  briefly  discuss  the  log-normal  model.  According  to 
Dyer  [36]  it  is  the  most  universally  applicable  distribution.  However, 


Litvinov  [39]  finds  it  to  be  non-applicable  to  rain  from  melting  snow, 
sleet  and  hail.  He  later  [40]  shews  explicitly  that  it  does  not  fit 
data  ftrom  melting  hail.  Incidently,  according  to  [40],  the  log-normal 


distribution  was  suggested  by  Levine  142].  Because  of  Litvinov's  con- 
clusions we  do  not  consider  the  log-normal  distribution  further. 
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For  the  P-8  distribution 


2A/y3 


(3.3-3) 


We  shall  define  the  probability  density  function  for  a drop  diameter 
by 

nfdfl) 

p(d0)«— . (3.3-1+) 

For  the  M-P  and  P-S  models 
-Ad. 

p(dQ)  « Ae  , M - P , (3.3-5) 

3 « -yd0 

P(d0}“2  d0*  U » P-S  . (3.3-6) 

< Since  p(dQ)  is  a function  of  rain  intensity  I through  A and  y, 
it  must  be  considered  as  a conditional  density  if  I is  random. 
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U.O  ORDINARY  SYSTEM  PERFORMANCE 

When  a linearly  polarized  radar  operates  in  rain,  using  a frequency 
above  a few  gigahertz,  there  may  be  undesirable  clutter  from  backscatter 
due  to  raindrops.  The  usual  way  of  reducing  this  clutter  is  to  resort 
to  circular  polarization.  If  the  system  transmits  and  receives  perfectly 
circular  polarization  in  all  directions,  and  if  raindrops  are  all  spherical, 
the  clutter  may  be  totally  cancelled  in  theory.  Neither  of  these  conditions 
is  satisfied  in  reality  and  the  purpose  of  this  section  is  to  determine 
approximate  bounds  on  the  practical  clutter  cancellation  achievable  by 
a system  designed  to  radiate  and  receive  approximately  circular  polar- 
ization. Practical  effects  of  raindrop  shape,  size  distribution,  and 
drop  orientation  are  determined  for  various  rain  intensities  from  0.5 
to  128  mm/hr. 

Two  situations  are  developed.  First,  the  raindrops  are  assumed 
to  be  perfectly  spherical  and  the  effect  of  system  polarization  tolerance 
is  found.  Second,  the  system  polarization  is  assumed  perfectly  circular 
and  the  limits  due  to  practical  rain  are  found.  In  the  latter  case, 
both  the  Marshall-Palmer  (M-P)  and  Polyakova-Shifrin  )P-S)  raindrop 
models  are  used  to  achieve  numerical  data. 

l+.l  System  Model 

The  system  model  is  that  shown  in  Figure  2.5-1,  and  described 
in  Section  2.5.  Th?  matrices  ft]  and  [r]  are  chosen  to  give  orthogonal 


space  polarizations 
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4.2  System  Clutter  Powers 
Ordinary  System 


Cancellation  ratio  CH  will  be  defined  as  the  ratio  of  clutter 
powers  *>ref/Por<j  where  Por^is  the  output  clutter  power  of  an  ordinary 
or  circularly  polarized  system.  P ^ *s  the  clutter  power  of  a "refer- 
ence" system  defined  as  having  a single  polarization  (say  X),  the  same 
pattern  and  transmitter  power  as  the  ordinary  system  and  operating  in 
the  same  rain  environment. 

Using  complex  notation,  PQrd  is  the  expected  value  of  |er|2  in 
Figure  2.5-1,  where  we  henceforth  drop  the  functional  dependence  of 
all  quantities  on  k since  we  deal  only  with  a typical  cell.  Now  er  = 
el  + e2*  where  ei  e2  are  Produced  by  the  sums  of  the  backscattered 
fields  from  the  N particles  in  cell  k.  Carrying  out  the  necessary 
algebra  leads  to  the  ordinary  system  power 


ord 


|A|*  2 { l |G(a  , 6 )|J 
i»l  1 


Kli  ^J2TB12i 


T2s 


22i 


(4.2-1) 


Here 


A 


KeT 

✓5  z2 


(4.2-2) 


Z is  the  radial  distance  to  the  center  of  cell  k,  K is  a constant  re- 
lated to  the  range  equation,  6 « 2*/A,  \ is  wavelength,  and  J * /-l. 

In  arriving  at  (4.2-1),  use  has  been  made  of  the  fact  that  s.^ 

* Also,  it  has  been  recognized  that  radial  position  of  raindrops 

is  unifromly  distributed  over  a range  cell,  which  is  large  relative 
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1 
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to  A,  and  position  is  independent  of  all  other  randan  quantities.  Hence 
the  expected  value  with  re3pect  to  radial  drop  position  has  already 
"been  taken. 

Reference  System 

The  reference  system  clutter  power  is  most  easily  obtained  by 
setting  s^2^  = 0 and  = 0 in  (l*.2-l)  and  adding  a factor  2 to  ac- 

count for  the  fact  that  the  two  systems  must  radiate  the  seme  total 
power 


U.3  Clutter  Cancellation  Ratio 

Clutter  cane  llation  ratio  CR  becomes  the  ratio  of  (U.2-3)  to 
(t.2-l).  The  expression  simplifies  by  observing  that  the  raindrop 
location  angles  and  6 may  be  considered  statistically  independent 
of  all  other  random  quantities  which  are  related  to  drop  scattering 
properties.  It  is  furthermore  reasonable  to  presume  the  same  statis- 
tical distribution  function  for  these  angles  applies  to  all  drops. 

With  these  considerations  we  have 

2*1  |0ll|2> 

CR  = — jj — . (k.3-1) 

I'll  -+J2TS12  - A?/’ 

Further  simplification  follows  the  assumption  that  N is  nonrandom. 

Such  a situation  is  true  if  +ho  rain  intensity  I is  taken  as  constant,. 
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Making  this  assumption,  and  noting  that  the  same  statistical  distri- 
bution applies  to  all  raindrops,  we  get 


CR  » 


2 M|au|2> 

«|.nlj2T  .^-T2  ,22|2  ) 


(U.3-2) 


Here  we  have  dropped  the  subscript  denoting  particle  i since  it  is  no 
longer  required.  We  consider  two  special  cases  of  (U.3-2). 


Ideal  Rain  Case 

For  spherical  raindrops  s^g  “ 0 and  s^  » Sgg.  Cancellation 
ratio  becomes 


CR 


2i2 


U - * I 


(fc.3-3) 


This  result  describe?  the  performance  of  a nonideal  system  oper- 
ating in  light,  drizzle-type,  rain  since  drops  become  more  nearly 
spherical  as  rain  intensity  decreases.  Figure  U. 3-1  is  a plot  of  CR 
from  ( 1* . 3—3 ) as  a function  of  |t|  with  thephase  of  T as  a parameter. 
To  obtain  30  dB  of  cancellation,  a system  must  not  have  more  than  0.2 
dB  of  amplitude  unbalance  in  the  two  linear  polarization  components 
for  zero  phase  unbalance.  Phase  unbalance  rapidly  decreases  the  al- 
lowable amplitude  tolerance. 

Ideal  System  Case 

It  can  he  shewn  that  the  functions  s _ are  all  real.  With  T = 1 


(^.3-2)  simplifies  to 


CANCELLATION 


B - -i!C1(lO”9)/6, 


(U. 3-7 ' 

(U.3-8) 
(^.  <-9> 


f 
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e.  is  the  dielectric  constant  of  air,  e is  the  dielectric  constant 
^ r 

of  rainwater  (relative  to  air),  and  dg  is  in  millimeters.  The  angles 
C and  iji  define  the  f all- path  of  the  raindrop  with  t,  being  the  azimuthal 
angle  of  the  vertical  plane  containing  the  trajectory  of  the  drop  while 
♦ is  the  fall  angle  measured  from  vertical. 

It  is  quite  difficult  to  solve  (k.3-5)  exactly.  Part  of  the  dif- 
ficulty arises  "because  C and  t|)  are  complicated  nonlinear  functions 
of  wind  components  and  4*  is  additionally  a function  of  dg  through  the 
drop  terminal  velocity.  In  general,  the  winds  should  be  treated  as 
random  quantities.  Other  difficulties  derive  from  the  complicated 
behavior  of  g^  and  g^.  Most  of  these  problems  can  be  avoided  if  we 
seek  to  calculate  bounds  on  performance,  rather  than  calculate  CR  ex- 
actly. 

For  the  lower  bound  on  CR  we  may  show  by  direct  calculation  that 
(g^  - gy)  is  positive  while  gy  is  negative  and  |gx  - gy|  _<  |gy|  for 
all  dg.  Thus,  as  a function  of  c,  (k.3-5)  is  minimum  when  sin  (c)  = 1. 
Presuming  this  to  be  true  for  the  moment,  the  resulting  expression 
is  minimum  as  a function  of  t|>  when  sin  (ip)  = 1.  We  define  the  final 
function  as  CR^^.  It  Is  given  by 


, 8 

”dn  s«Sl-g)5) 


More  careful  study  of  the  behavior  of  x,  and  iji  shows  that  (U.3-10)  re- 


sults if  we  make  the  simple  assumptions  that  all  winds  are  constant 
(nonrandom)  and  that  a strong  horizontal  wind  blows. 
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For  the  upper  bound,  ( U. 3—5)  is  maximum,  as  a function  of  for 

2 

sin  (?)  ■ 0.  The  resulting  expression  becomes  maximum  as  cos(^)  approaches 
zero.  Again  assuming  nonrandom  winds  with  a strong  horizontal  compo- 
nent V , the  upper  bound,  defined  as  CP  , becomes 
r * **  * max  * 


CR  » 
max 


2 VrU  B(gy2} 

E{(gK-gy)?(vl-Vv)S 


( U. 3-11) 


Here  V is  the  vertical  wind  speed  and  is  the  drop  terminal  velocity, 
which  is  approximately  given  by 

Vm  = 10.105[1  - e~*°/2)  (U.3-12) 

if  dQ  is  in  millimeters. 

A final  situation  supposes  all  winds  are  zero.  The  (22)  reduces  to 


a(«r2> 

E{(ex-gy)2) 


(no  winds) 


(U. 3-13) 


U.U  Ideal  System  Numerical  Results 

Calculation  of  (U.3-10),  ( U . 3—11 ) and  ( U . 3—13 ) is  still  not  a 

simple  task  due  to  the  form  of  the  expressions  for  and  g . However, 

2 2 2 

direct  calculation  of  the  functions  g , g and  (g  -g  ) has  shown: 
l)  that  sensitivity  to  the  parameter  is  small  and,  2)  they  may  he 
approximated  closely  with  simple  functions  which  lead  t,o  closed-form 
solutions  for  the  various  expectations.  Thus,  approximations  were 
developed  and  expectations  were  taken  assuming  ra  ndrops  distributed 
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according  to  both  the  M-P  distribution  and  the  P-£  distribution. 
Marshall-Palmer  Distribution 

The  applicable  distribution  probability  density  is  (3.3-5)  where 
A is  given  by 

A ■ it.ll”0*21  (U.fc-l) 

where  I is  rain  intensity  (mm/hr).  With  this  density  function  (U.3-1C), 
(It. 3-11)  and  ( it . 3—13)  evaluate  to 


CR 


min 


122.77  A[1.U05  - 1.656A  + A3! 
155-63  + A2 


123.96  [O.J53  t..\gl 
^ + (155.63X,2/A2) 


(U.U-2) 


(U.U-3) 


CR 


no 


= 123.96  A2[0.7;3  ♦ A2] 
155.63  + A2 


(lt.U-U) 


Assuming  V « 0.  Here 

w it  + 6 it  

1 = " [1  + (0.5/X)]9  + [1  + Cl/A)]9  [1  + (1.5/A))9 


[1  + (2/A)]9 


(it.lt-5) 


[1  + (0. 5/> ) JXl  [1  + (l/A)]J"L  [1  + (1.5/A)] 


11 


,11 


(1  + (2/A)]11 


(it.  It— 6 ) 
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Figure  4,4-1  illustrates  plots  of  (4.4-2),  (4.4-3)  and  (4.4-4). 
Polyakova- Sh.ifr  in  Distribution 

Here  the  probability  density  function  of  drop  diameters  is  given 
by  (3.3-6)  and  y is  given  in  'Table  3.2-1.  Calculation  now  produces 


CR 

max 


Vr  ] 77.21  (1.21  + y2) 

i10,105/  + (228.467X2/y2) 


(4.4-7) 


CR 


min 


76.469  y (2.76  - 2.662  y + v3) 
228.467  + Y2 


(4.4-8) 


CR  (no  wind)  = » ft 

228.467  + y2 


(4.4-9) 


where 


X!  " 1 - 


[1  + (0.5/Y)] 


11 


[1  + (1/y) ] 


-11 


[l  + (l. 5/y) ] 


li 


+ rr  (4.4-10) 

11  + (2/y)]11 


4 + 6 4 

X?  = 1 " (1  + (0,5/y)]13  [1  + (1/y)]13  [1  + ( 1 . 5/y ) ]13 


+ - 


1 

13 


h + (?/y)1 


(4.4-11 ) 
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The  expresaions  (i».U-7)  through  (U.U-9)  are  plotted  in  Figures  U.lt-g 
through  for  rain  from  melting  snow,  sleet  and  hail. 
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5.0  PERFORMANCE  OF  CANCELLATION  SYSTEM  WITH  UNIFORM  RAIN  CHARACTERISTICS 

5.1  Rain  Characteristics 

As  in  section  3 we  consider  rain  particles  to  be  nonspherical  in 
this  section.  However*  unlike  chapter  3,  we  shall  here  assume  particles 
all  have  uniform  shape,  size,  and  orientation. 


5.2  System  Equations 

The  cancellation  system  is  shown  in  block  iiagiam  form  in  Figure 
5.2-1.  At  the  left  of  the  diagram,  the  transmitter  signed  is  split 
in  a power  sense.  Half  of  the  signal  power  is  directed  to  the  antenna 
which  generates  the  X component  of  the  radiation  field,  dhe  other  half 
of  the  signal  power  is  directed  through  the  lower  channel  and  eventu- 
ally is  used  to  generate  the  Y component  of  the  radiation  field.  How- 
ever, this  second  portion  of  the  signal  is  modified  in  magnitude  and 
phase  by  the  element  T.  It  is  this  element  T which  represents  the  basic 
system  imperfections  which  prevent  the  transmission  of  a truly  circu- 
larly polarized  field.  Finally  the  signal  in  the  Y channel  is  shifted 
in  phase  by  + ir/2  with  the  +v/2  corresponding  to  polarization  of  a lef* 
hand  sense  and  the  -n/2  producing  polarization  of  a right  hand  sense. 

The  quantities  ■>  and  e * * represent  the  electric  field*  produced 

* y 

by  the  X and  Y antennas  respectively.  By  using  equation  (2.U-1),  one 
may  express  these  fields  as  indicated  in  equation  (5.2-1). 


e 

/2  Z 


(5.2-la) 


♦1TK,emG„ 
- ITT 


-J80Z 


(5.?-lb) 


f 


i- 

f 


t, 

I 


1 

r 

\ 
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If  the  transmitter  antennas  are  not  perfectly  orthogonal  then  a 

cross  coupling  of  the  electric  field  components  occurs  as  indicated 

by  the  t^^  factors  in  Figure  5.2-1.  This  cross  coupling  is  expressed 

mathematically  in  equation  (2.5-1),  and  the  t . , terms  are  defined  in 

equation  (2.5-2).  Thus,  the  true  orthogonal  components,  e+  and  e , 

x y 

may  be  expressed  as 


- *11 V * ‘l2ey 


and 


( 5 • 2~2a ) 


+ ^ + + 

e - t.,e  - + t„„e  ■ 
y 21  x 22  y 


(5.2-2b) 


The  relationship  between  the  incident  field  components,  and 
e+  , and  the  reflected  or  scattered  field  components  is  given  by  equation 

y 

(2.3-3)  with  the  terms  defined  in  equation  (2.3-6).  These  equations 
together  with  equation  (2.U-2)  which  indicates  the  propagational  mag- 
nitude and  phase  effects  may  be  used  to  express  the  scattered  field 
components  at  a distance  Z'  from  the  particle  as 


- ki<  +si2%)  -jv' 


(5-2-3a) 


e - 

y 


(S2lex  * S22ey^  "J80Z 


(5. 2- 3b) 


As  in  the  case  of  the  transmitter  antennas,  the  axes  of  the  re- 
ceiver antennas  may  not  be  perpendicular  to  one  another  so  that  cross 
coupling  of  the  field  components  occurs.  The  matheraot ica]  basis  for 


this  effect  is  given  in  equation  (2.5-^),  and  the  representation  used 


1 
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in  this  work  is  shown  in  equation  (5.2-4). 


V 3 ruex  + r12ey 


(5.2-4a) 


V=r21ex  +r22ey 


(5.2-Uh) 


Receiver  signals  e^  and  e^  may  be  determined  by  application  of 
equation  (2.4-3). 


el  = h2Gr'x'  ’ K2tV[rllex  * r12ey1 


K2°r  -ll(3llex  * * r!2(S2iex  't  822*v0 


z' 


K2Gr  [rll(Sll  (*11 V * ^v)*  ^(^lV  * t22ey^)) 


+ rl2  ^S21 


/ + + \ / + + \ , ) 
VllV  + ^y-)  + 8»\Vx'  +t22fty7  >J  e 


W r , 

[rll(Sll  \ tll°T  iJ,rt12GT/  + S12  iJTt22GT,  > 


r12  ^S21  {^l0?  lJTt12GT  ) + S22  (t21GT  iJTt22GT  ) ll 


-je0(z*z') 


(5. 2- 5a) 
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2 ^2  Z«Z' 


r21(Sll  <*U°?  iJTt'lS°T  )>  * S12  )> 


+ r 


22 


(821<V?  * S22  (*21$  iJTt22GT>)]e"JB° 


(z+z') 


(5. 2- 5b) 


For  the  case  of  a monostatic  radar  the  following  relationships 
are  valid. 


Z = 7/ 


R - T 

X _X 
eT  GR 

Y Y 
GT  “ °R 


= G 


5 G' 


The  assumption  that  the  antenna  axes  are  orthogonal  results  in  a sig- 
nificant simplification  in  equations  (5.2-5a)  and  (5.2-5b).  This  sim- 
plified form  is  given  in  equations  (5.2-6a)  and  (5.2-6b). 


el  - 


„ w 


f2  if 


[o\x  +JTGYS121  e 


-J2S0Z 


(5.2-6a) 


e2  “ 


+JRK1K2eTGi 


/?  z‘ 


[c 


gxs21  .jtgys22 


]• 


-J2ecz 


(5.2-6b) 


Use  of  the  definition  in  equation  (P.U-.U)  along  with  the  assumption 

X Y 

that  the  X and  Y radiation  characteristics,  G and  G , are  Identical 
results  in  a further  simplification  for  the  terms  representing  the 
signals  received  an  a result  of  scattering  by  one  particle. 


f 
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Ke^fo,  6) 

" /2Z2 


[sii 


-J2B0Z 


, and 


(5.2-7a) 


♦JTKeTG2(a,  6) 
«/2  Z2 


[S21  +JTS22]  e'tJ260Z 


(5.2-Tb) 


where 

GX  = 0Y  « G(a,  0)  . (5.2-8) 

In  the  process  of  evaluating  actual  radar  system  performance  in 
the  presence  of  rain,  it  will  be  necessary  to  consider  multiparticle 
scattering.  As  pointed  out  in  section  2.U,  the  multiparticle  effect 
is  evidenced  by  a summation  of  individual  particle  returns  such  as  those 
of  equation  (5.2-7).  Expressions  for  e^  and  e2  for  multiparticle  scat- 
tering with  the  same  assumptions  as  used  in  equations  (5.2-7a)  and 
(5. 2- 7b)  are  given  in  equations  (5.2-9a)  and  (5.2-9b). 


e1(k) 


KeTe 


-J£ezk  N, 


6i)[slli(k)  iJ«12l(*)] 


-J2VZi 


(5.2-9a) 


e2(k) 


lJTKeTe 


-JS8Zk 


i 02(«It  B,)  [s. 

i“l 


21i(k)  +jTS?2i(k) 


•J2BoAZi 


( 5 . 2-9b ) 


Before  it  will  be  possible  to  continue  the  analysis  of  the  can- 
cellation system,  one  must  define  the  G element  in  the  system  block 
diagram.  The  original  definition  of  G was  provided  by  Rice  and  Peebles, 


and  it  is  provided  brlov. 
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, e At) 

G(k)  " TaTTiT  X,  ^*7  (5.2-10) 

fc*k-M 

The  value  of  G at  the  time  the  system  is  processing  the  returns  from 
th 

the  k— 1 - range  cell  is  determined  by  the  negative  of  the  average  value 
of  the  ratio  e^/e,,.  This  ratio  is  evaluated  in  the  M cells  before  cell 
k,  cell  k itself,  and  the  M cells  after  cell  k. 

Using  the  definitions  of  e^  and  provided  by  equations  (5.2-9a) 
and  (5.2-9b),  one  obtains  the  following  expression  for  G. 


N. 


£ G-(«i,8i)[Slli(i)  +JTS12i(i)] 
= +JT(2M+1)  Z “n“ 


-J20OAZ. 


(5. 


£=k  _M  ^ o ~J28_dZ 

£ G («!»  8.)[S2i.(Z)  +jTS?2i(H)]  e ° 1 

i=l 


If  all  the  particles  within  a given  range  cell  have  the  same  sice, 
shape,  and  orientation  then 


G(k) 


k+M 

+JT(2M  +1 ) Z 
£=k-M 


[Si;l(£)  +JTC12ff.)] 

TOrr^CTirr 


(5.. '-12) 


G may  be  expressed  by  a yet  simpler  form  by  requiring  cell-to-cell 
uniformity  in  all  particle  parameters.  This  assumption  gives 


G(k)  = 


-1 

+jf(2M+l) 


( 2M+1 ) 


t5n  iJTSig1  -1  lsn 

Is?1  +JTS??]  ' +JT  (s?1  +JTS,0] 


(5.2  -l'<) 
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All  of  the  equations  developed  thus  far  in  this  section  pertain 


to  the  situation  in  which  only  clutter  is  present  in  the  cells  being 


observed.  In  the  event  that  a target  appears  within  one  of  the  cells 


of  interest,  one  obtains  the  following  expressions  for  e^,  e^,  and  G. 


e1(k) 


e£(k)  + e1(k) 


(5.2-lU) 


e2(k) 


e^(k)  + e^(k) 


(5.2-15) 


The  quantities  e^(k)  and  e^Ck)  are  the  clutter  contributions  to  the 


total  backscattered  signal  and  will  in  general  be  expressed  in  the  form 


indicated  in  equations  (5.2-9a)  and  (5.2-9b).  Target  contributions 

T T/ 

to  the  backscattered  signal  are  indicated  by  e^(k)  and  e^(k)  and  these 
terms  find  expression  in  equations  (2.5-17a)  and  (2.5-17b).  Application 


of  equations  (5.2-lU)  and  (5.2-15)  to  general  definition  of  G in  equation 


(5.2-10)  gives 


G(k)  = 


k+M  Vn  ei'-' 

-1  ei  -1  ^ ~7tt  C1 

(2M+l)  Z-/  elu)  ~ (2M+1)  l- k-M  *2  ’ + 

i=k-M  Ifr 


eAD 


-K-j  eq(k) 
2 * + ^TTkl 


(5-2-16) 


Finally,  using  the  appropriate  combination  of  expressions  for 


e^v  e^,  and  G,  one  may  express  *,he  receiver  output,  e^ , in  general. 


er(k)  =•  e1(k^  + G(k)e2(k) 


(5.2-17) 
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5.3  Evaluation  of  Performance 

In  this  section,  the  output  of  the  cancellation  system  as  indicated 
generally  in  equation  (5.2-17)  will  be  evaluated  for  several  different 
situations.  These  situations  are  itemized  below. 

1.  Only  clutter  is  present  in  the  radar  system’s  field  of  view. 

2.  Both  a target  and  clutter  are  in  the  system's  field  of  view, 
and  target  occurs  in  the  range  cell  presently  bei,.g  examined. 

3.  Both  a target  and  clutter  are  in  the  system's  field  of  view. 

The  target  occurs  in  a range  cell  whose  returns  are  used  in 
the  evaluation  of  he  present  G function,  but  that  cell  is 
not  the  cell  presently  being  examined. 

For  the  first  case,  the  needed  expressions  for  e^,  e^,  and  G are 
provided  by  equations  (5.2-9a),  (5.2-9b),  and  (5.2-13).  The  receiver 
output  signal  is 

e (k)  = e,(k)  + G(k)  e„(k).  (5.3-1) 

T X e 

Equation  (5.2-13)  reflects  the  assumption  of  cell-to-cell  uniformity, 
and  when  this  concept  is  applied  to  the  definitions  of  and  e,j  a sim- 
plification in  form  results. 


e,  (k)  = 


KeTe 


-J28ozk 


Jz  z‘ 


[SU  :JTS12 


\ 

)£ 

i=l 


(«.  , Bi  )e 


-J2Vzi 


( 5. 3-2a) 


,(k)  » 


^jTKeT 


-^2eozk 


[S?1  +.1T3??1^  G‘(&i 
1 = 1 


B )r 


(5. 3-.?b) 
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G(k) 


-1  ^Sll 

+JT  fs21  +JTS223 


(5.3-3) 


One  may  now  evaluate  the  system  output  under  the  assumption  of  cell- 
to-cell  uniformity  by  substituting  the  above  expressions  for  e^,  eg. 
and  G into  the  system  output  equation,  (5.3-1). 


N. 


er(k) 


Ke„e  " " **  ~J2enAZ, 

~-~2 [SV  +JTS12]  £ G Me  01 

S2Z?  ^ i»l  1 1 


^ /-I  \ tSll  tJTSl2] 

\-JT/  ^S21  ^JTS22^  ST Z ‘ 


[S21  i'1TS22] 


y -J20  AZ 

•£  G <v  Ve  = o 

i=l 


(5.3-1*) 


In  the  second  case,  the  appropriate  expressions  for  e^,  e,. , and 
G come  from  (5.2-1**),  (5.2-15)  and  (5.2-16). 


e (k)  = e1 (k)  + G(k)e0(k) 
r i - 

= e^k)  + ^M+1j 

ay  k 

- (*Kfc>  * Mr) 

tyk 


k+M 


l i=k-M 


ex(.e) 


,(k ) 


k+M 

Ei: 

|t=k-M 


e^U) 

e“CI7 


ex(k) 

+ eJkT 


e2(k) 


(5.3-5) 


(5.3-6) 
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e1(k) 


(S  1JTS  ) 
—±Lz i£_  e (v) 

Ts^+iTs  ) VK; 


(5-3-7) 


(5.3-8) 

What  seams  to  be  a new  parameter,  G^k),  appears  in  the  last  step  of 
equation  (5.3-8).  Actually,  comparison  of  G^(k)  with  the  G(k)  of  equation 
(5.2-13)  leads  one  to  conclude  that  the  two  parameters  are  the  same. 

This  latter  observation  simplifies  evaluating  equation  (5.3-8)  since 
the  first  part  of  the  equation, 

e^(k)  + G^kje^k)  , 

should  now  be  recognized  as  being  the  same  as  equation  (5.3-b).  Thus, 


er(k)  = 


[e^(k)  + G1(k)e^(k)l 


(5-3-9) 


The  statement  defining  the  third  case  specifies  that  while  the 
target  is  not  in  the  cell  presently  being  examined  it  is  located  in 
one  of  the  cells  being  used  to  evaluate  0.  It  is  assumed  for  this 
analysis  that  the  target  is  in  cell  m,  and  thus  e^(k)  and  &g(k)  are 
given  by  equations  (5.3-2a)  and  (5.3-2b)  respectively  while 


G(k) 


(2M+1) 


k+M 

£ 

t=k-M 

Ifm 


e U) 


el(m) 


e'2 ( i)  * e0(m) 


(5.3-10) 


This  expression  for  G along  with  the  specified  expressions  for  e^  and 
e0  may  be  used  to  obtain  the  following  expression  for  e . 


79 


er(k)  * e1(k)  + G(k)e2(k)  ■ ^ 
e^w)  e^(m)  + e^(m) 


e“(k) 


e2(k) 


e^un;  C,  % ^ T,  « 

2 e2(m)  + e£(m) 


(5.3-11) 


(5.3-12) 


Equation  (5.3-11)  will  be  further  analyzed  for  two  special  situations. 
The  first  situation,  which  might  be  referred  to  as  the  high-noise  en- 
vironment, corresponds  to  the  mathematical  statements 


e^(m)  >>  e^(m) 


and 


(5.3-13a) 


e2(m)  >> 


(5.3~13b) 


When  equations  (5.3-13a)  and  (5.3-13b)  are  applied  to  equation  (5.3-12), 
the  following  approximations  are  valid. 


e^(m) 

e2(m) 


(5.3-10 


(5.3-15) 


The  following  results  are  useful  in  working  with  equation  (5.3-15). 


ej(m) 

e^m) 


/ 1 \ ^Sll  -JTS12^ 

\+jt;  r s21  ?jts?2] 


(5.3-16) 


ry  ihhjt^vi;  « 
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j— i "fTvpr'/j  wyr~’r’r,^-ryy  ^t’Y  ^^*7^ 


(&)  ••“ ' 


”J260Zk  N. 

TKeTe  11  K 


*4 


^ 2 -J28 AZ 

fSn  +JTS12]  2^  G2(a.,  Si)e  0 1 
1=1 


ej<k) 


(5.3-17) 


The  result  of  applying  equation  (5.3-17)  to  equation  (5.3-15)  is  given 
in  equation  (5*3-18). 


.(*)  “ (2M+1)  fei(k)  " e^(k))  = 0 


(5.3-18) 


A second  interesting  situation  arises  for  equation  (5.3-11)  when 


T C 

e1(m)  » e^(m)  and 


(5.3-19a) 


e2(m)  >>  e2(m) 


(5.  3-  19b) 


These  assumptions  lead  to  what  might  be  called  the  large  signal  environ- 
ment, and  the  output  signal,  e , is  determined  to  be 

r 


er( 


" (am) 


c /e^(m)\ 

h(k)'w^v 


e?(k) 


(5.3-PO) 


5.**  Comparison  of  Performance 


In  section  5*3,  expressions  were  developed  for  the  output  of  the 
cancellation  system  under  various  conditions,  and  in  this  section,  the 
results  obtained  for  the  cancellation  Bystem  shall  be  compared  with 
the  appropriate  results  for  the  reference  system.  The  basis  of  com- 
parison, in  most  cases,  shall  be  the  ratio  CF  which  is  defined  below. 


«j 


J 
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CP  . output  power  of  the  reference  system  < , 

output  pover  of  the  cancellation  system  0.4-1 

The  method  of  power  calculation  requires  some  discussion.  Figure 
(5.k-l)  shows  a potential,  E,  applied  across  the  terminals  of  a general 
impedance,  Z.  From  the  theory  of  netovrk  analysis  [h3l,  one  recalls 
that  when  E and  Z are  represented  in  phasor  notation  the  complex  pover 
delivered  to  or  absorbed  by  Z i s 


In  the  present  analysis,  E represents  the  output  signal  of  the  radar 
receiver,  and  Z corresponds  to  load  impedance  of  the  network  driven 
by  the  receiver  output.  Now  apply  the  following  definitions. 


er^REF  = 


ZREF 

PREF 

er(k) 

ZCAN 


CAN 


the 


the 


the 


PCAN  = the 

cR-^s; 

PCAN 


output  signal  from  reference  receiver, 
load  impedance  of  the  reference  receiver, 
output  power  of  the  reference  receiver, 
output  signal  from  cancellation  receiver, 
load  impedauce  of  the  cancellation  receiver, 
output  power  of  the  cancellation  receiver. 


(5.U-3) 


^REF ^ 


l*r *k)CAN ' 


ZCAN 

ZREF 


(5.h-U) 
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If  ZCAH  * W then 


CR 


ler^REF^ 

ler*k^CAlJ 


(5.U-5) 


It  is  the  form  of  CR  indicated  in  equation  (5.^-5)  that  will  he  evalu- 
ated in  the  following  material. 

Equation  (5«3-^)  shows  that  the  output  of  the  cancellation  system 
is  zero  when  only  clutter  appears  in  the  field  of  view  of  the  radar, 
and  thus 


l'r(lIW2  ’ 0 


(5.k-6) 


The  output  of  the  reference  system  given  the  same  constraints  imposed 
on  equation  (5.3-*0  is 

N. 


Ke  s -J28-Z,  * 2 . 

M-W  ■ -pr1  • i'  i 

k 


-3 2B0AZi 


(5  U-T) 


(5.U-8) 


Thus, 

| «r )rEF I ^ " °2  **  0 » and 


(5.M) 


the  cancellation  ratio  may  be  shown  to  be 

«•#-- . 


(5.*»-l0) 


The  significance  of  equation  (5«h-10)  is  that  it  states  that  the 
cancellati • n system  offers  Infinite  improvement  in  clutter  cancellation 
when  compared  to  the  reference  system. 
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The  radar  system  output  under  the  conditions  outlined  in  case 
two  of  section  5.3  contains  hoth  a desired  component , referred  to  as 
the  signal,  and  an  undesired  component,  called  noise.  When  the  system 
output  ccntains  hoth  desired  and  undesired  components,  one  finds  that 
an  excellent  criterion  for  Judging  system  performance  is  a ratio  which 
indicates  the  relative  strength  of  the  signal  to  the  noise.  This  cri- 
terion is  referred  to  as  a signal-to-noise  ratio  [M],  and  the  ratio 
is  defined  in  equation  (5. ^-11). 


SNR  = 


output  si< 


output  noise  power 


(5.1+-11) 


The  reference  system  has  an  output 


e (k)  * e^(k)  + eT(k)  (5. b-12) 

r 

C T 

where  e (k)  is  the  noise  or  clutter  component,  end  e (k)  is  the  signs] 

or  target  component.  This  leads  to  the  signal-to-noise  ratio  for  the 

reference  system  given  in  equation  (5- *+-13). 


It  is  evident  from  equation  (5.3-9)  that  the  cancellation  system 
output  contains  no  noise  or  clutter  component  so  that  the  signal-to- 
noise  ratio  for  the  cancellation  system  is  infinite. 

When  the  operating  conditions  satis fy  the  specifications  of  the 
third  situation  in  section  5>3,  one  wishes  to  minimize  the  system  out- 
put since  any  output  under  these  conditions  amounts  to  noise.  The 
cancellation  ratio  will  again  be  used  a the  criterion  of  performance, 
and  the  following  expression  . &■*<?  used  in  evaluating  the  cancellation 
rati 


l 
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er^REF 


■Vu  2,  . ^260AZi 

2-G  (“i*  si)e 


i*l 


ASuB 


(5.U-1U) 


KeT 

A * -5-  e 

4 


(5-*»-15) 


N. 


* -J20  AZ 

B = X)  G2(ai»  Bj)© 


(5.U-16) 


i«l 


AB 1 2 8^ 


( 5 . **-17 ) 


The  output  of  the  cancellation  system  in  a high  noise  environment 
as  given  hy  equation  (5*3-18)  is  zero,  and  thus,  the  output  noise  power 
is  zero.  Thus, 


CR  - 


(5.^18) 


If  the  large  signal  environemnt  assumptions  are  used,  one  finds 
that  the  output  of  the  cancellation  system  is  given  by  equation  (5.3-20). 
The  notation  may  be  simplified  if  the  following  definitions  are  applied. 


e^(k)  « 

ABC//2 

(5.  **-19) 

“ 

+JTADB//F 

( 5.  **-20) 

C“S11 

-JTS12 

(5.  **-21) 

D,  S2i 

♦jts22 

(5.  '*-22) 

F * <*T(m)/e2(m) 

(5. *'-r3) 
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The  cancellation  system  output  in  the  above  notation  is 


AB 


“r<k>C»  ' Jz  (2M+1) 


[C  +JPTD] 


^(k)CAlJ  * ',o  MCI2  +2Im(C*FTD)  + |FTD|2] 


2(2M+l)‘ 


(5.U-2M 

(5.U-25) 


The  large  number  of  parameters  in  this  equation  makes  it  difficult  to 
see  the  nature  of  overall  variations.  The  following  assumptions  are 
restrictive  in  nature  but  do  aid  in  clarifying  the  system  effects. 

T « 1.0  (5. fc-26) 

F * * cob(«)  + J sin(0)  (5. ^-27) 


One  further  point  bears  mentioning  before  the  analysis  is  continued. 
Since  the  output  of  both  the  reference  system  and  the  cancellation  system 
is  randan  in  nature,  theoretical  correctness  demands  that  the  expres- 
sion used  for  output  power  be  a statistical  expectation.  Thus  equation 
(5.^-5)  becomes 


CR  ■ 


«Ur(k)CAll|2) 


where  E{  • } denotes  a statistical  expectation  [ 15 ] . 


( 5. k-28) 


E{|er(k)REF|2>  = E{  |ab|2}  * (5.«*-29) 

E{  I e (k)p.J2}  = E{  • [|c|2  +2Im(C*DeJ0)  + |d|2]  (5.U-30) 

r CAW  2(2M+1)2 


■ . MWI 


mmrnimmmmtQim 
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2S2^(2M+l)2 

[ | C | 2 +2£n(C*Deje)  + |d|2] 


(5*l»-3l) 


(5.1*-  32) 


An  interesting  result  is  obtained  if  the  cancellation  ratio  is 
evaluated  with  the  output  of  the  ordinary  system  replacing  the  output 
of  the  reference  system  and  if  0 equals  180  degrees.  Under  the  stated 
conditions 


CR  - (2M+1)2 


(5.U-33) 
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6.0  PERFORMANCE  OF  THE  CANCELLATION  SYSTEM  WITH  A REALISTIC  RAIN  MODEL 

6.1  The  Rain  Model 

There  are  basically  two  major  areas  involved  in  specifying  a rain 
model.  The  first  area  is  concerned  with  describing  the  individual  rain 
particle,  and  this  aspect  of  rain  modeling  has  been  delt  with  in  chapter 
3 where  a major  emphasis  was  on  defining  particle  size  and  shape. 
Particle  size  and  shape  affect  the  radar  system  equations  through  the 
S matrix  elements  defined  in  equation  (2.3-6),  and  nothing  in  this 
chapter  will  be  done  to  alter  the  nature  of  the  influence  that  particle 
size  and  shape  have  on  the  system  equations. 

The  second  aspect  of  rain  modeling  is  concerned  with  the  nature 
of  large  scale  rain  systems.  Here  individual  particle  characteristics 
are  lost  and  only  overall  effects  are  noted.  In  chapter  5 it  was  as- 
sumed that  the  rain  system  was  characterized  by  overall  uniformity, 
but  the  present  goal  is  to  establish  a less  restrictive  and  more  re- 
alistic model  for  the  rain  system.  In  pursuit  of  the  stated  goal, 
the  following  assumptions  will  be  applied. 

1.  The  rain  particles  exhibit  a uniform  spatial  distribution 
over  the  range  cells. 

2.  Particle  3ize  and  orientation  will  be  treated  as  random  vari- 
ables, and  it  is  assumed  that  they  are  statistically  indepen- 
dent of  the  particle  position  parameters.  The  distribution 
of  particle  size  presented  in  chapter  3 will  be  used  here. 

3.  The  number  of  narticles  occuring  in  a range  cell  is  strictly 
random,  and  the  distribution  function  for  the  number  of  par- 
ticles will  be  determined  in  later  work.  It  is  assumed  that 


C/A 
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the  masher  of  particles  is  statistically  independent  of  par- 
ticle size,  shape,  orientation  and  location. 

U.  Although  several  parameters  relating  to  the  rain  are  random 
in  nature,  it  is  assumed  that  the  statistical  distributions 
describing  the  various  parameters  are  invariant  in  time  and 
space. 

5.  It  is  assumed  that  the  clutter  returns  from  one  range  cell 
are  statistically  independent  of  the  clutter  returns  from  any 
other  range  cell. 

6.2  System  Equations 

The  physical  system  bieng  analyzed  here  is  the  same  as  the  one 
presented  in  Figure  5.2-1,  and  so  the  definitions  for  e^k),  e2(k), 
and  G(k)  given  in  equations  (5-2-9a),  (5.2-9b),  and  (5.2-10)  respectively 
are  valid  here.  However,  for  several  reasons  including  both  ease  of 
notation  and  statistical  considerations  the  following  forms  will  be 
used. 


e1(k)  « 


(X_ ( k ) + JY  (k) 
V2  ' X 1 


(6.2-1) 


e2(k) 

G(k) 


^2(k)  + jY2(k) 


(-1)  kfM  ’1<‘>  /-l\ 

Tfc+M 

V 

«i(A)  e^(k) 

\2m+i/  t-£.M  ^rn  \2m+i; 

l»k-M 

-itY 

^UT  + ^TkT 

C'(k)  + 


e^(k) 

e2(k) 


(6.2-2) 


(6.2-3) 


The  quantities  X^,  Y , X^,  and  Y2  are  determined  by  the  following. 

jYn(lO 


Slli^  -^TS12i^  = ali^k^  e 


S2ii(k)  +JTS22i(k)  = agi(k)  e 


jY2i(k) 


(6.2-4) 


(6.2-5) 


20QAZi  >=  Q± 

"*  , 

X^k)  = X!  G2(ait  0i)ali(k)  cos^y^k) 

\ 

Y-Jk)  = 53  G2(oi,  0i)a11(k)  sin^rli(k) 

\ o 

X2(k)  = G‘L(ai,  0i)a2i(k)  cos  ^2i(k) 

Nk 

Yg(k)  * ^ G2(a^ , 0i)a2^(k)  sin(y2i(k) 


\) 


- 9j 


-»4) 


(6.2-6) 


(6.2-7) 


(6.2-8) 


(6.2-9) 


(6.2-10) 


Equations  (6.2-1)  through  (6.2-10)  may  be  used  to  express  the 
system  output  as  well  as  the  magnitude  squared  of  the  system  output. 


er(k)  * e1(k)  + G(k)  • e2(k)  • e 


/ 2M\  , . 0 (k)  ‘ e2(k) 

\2M+1/  elU;  { 2M+1) 


1^  + (2M+1) 


e (k) 

°'(k)  + 


e2(k) 


(6.2-11) 


|er(k)|  = er(k)  • e*(k)  = (2^1)  Je1(k)  | 2 - - — — Re^kje^k^  (k)) 

' ' ( 2M+ 1 ) 


^'(k)!2  |e2(k) 
( 2M+1 ) 2 


(6.2-12) 
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};L(k)|2  * [x2(k)  + Y2(k)] 


(6.2-13) 


>(k)|2  - LT.L2.21A1  ■ [X2(k)  + Y2(k)] 


(6.2-lM 


e^k)  e2(k) 


(k)  X2(k)  + Y1(k)  Y2(k)J 


+ J^k)  Y2(k)  - X2(k)  Yx(k)) 


(6.2-15) 


^ [X1(£)  + JY1(£)] 

{k)  " £^1M  ^jrTngTT  +jY2nT! 

£*k 

+ ^ Re(T)(xg(£)Y1(t)  - X1(£)YgU))  - In(T)  (x^t  )Xp(£ ) ♦ Y^pYgU)) 

i=k-M 
£*k 


k+M 

* E 


|T|2  (: 


X2(£)  + Y2(£)J 


£=k-M 

£j<k 


(^k+M 

E 

£=k-M 
v £?<k 

('  k+M 

£ 

£=k-M 


M2|x 

+ 

oi 

C\J  C\J 

Y2(£)j 

(6.2-16) 

Re(T) 

(xg(£)Y1(£)  - 

- X1(£.)Y 

2(‘>)  + 

Im(T)  (xg(£)Y1(£) 

- X1(£)Y?(£^ 

|t|2  | 

X2(£) 

Y2(£)J 

Re(T) 

(xi(£)X2(£) 

► Yi( £)Y 

2<l))  + 

Im(T)  (x2(£)Y1(£) 

- XjUjYgfOj 

I T i 2 

22U)  * 

Yg(£)' 

(6.2-17) 

The  relationship  given  in  equation  (6,2-18)  will  be  useful  in  later 


work  with  equation  (6.2-17). 
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(£  •■)’ 


N 


i»l 


N 


•?  * 


E E 

i-l  J-l 

i*J 


SiaJ 


(6.2-18) 


As  a result  of  the  random  nature  of  the  rain  as  outlined  in  the 
preceeding  section,  much  of  the  evaluation  done  in  this  chapter  will 
be  done  in  a statistical  sense.  This  means  that  power  will  be  defined 
as  E(|xj  } where  the  operator  E{  ♦ } is  the  expectation  operator  presented 
by  Thomas  [ 1*5 ] • Thus  the  output  noise  power  of  the  cancellation  system  is 


E(|er(k >|*} 


E{  |e^(k) \)  - 

11  1 (2M+1)2 


UM  — E(Re  <e*(k)e2(k)G'(k)^  } 

+ (2^1)  E(|G'(k)|2  Ie2(k)i2}  (6.2-19) 

Now  consider,  in  order,  the  three  expectation  terms  appearing  on 
the  right  hand  side  of  the  equality  in  equation  (6.2-19). 


E{|ei(k)|P}  * [X2(k)  + Y2(k)]| 


(6.2-20) 


Examination  of  equation  (5. 1+-15)  reveals  that  A is  not  random  in  nature 
so  that 


E{|ei(k)|2)  » [E{X^(k)>  + E{Y2(k)}]  . 


(6.2-21) 


From  appendix  A, 


E{X2(k) ) = E(Nk>  e{°  | E{a2}  , and 


(6.2-22) 


E(Y2(k ) } = E(Nk)  E{°  E(a2] 


(6.2-23) 
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Appendix  B shove  that 
K2 

E{NkJ  - KjVj^Ig  , (6.2-21*) 

and  appendix  C indicates  that 

EfG^o.B)}  « 0.70808  . (6.2-25) 

2 

Consider  the  remaining  unknown  factor  Eto^}.  From  equation  (6.2-U), 

“l  * lSll  ±JTS12 * ‘ (6.2-26) 

Thus 

“l  * Sll+  lT!2  sx2  *2In^TSllSl^  » ***  (6.2-27) 

E{a*}  • EtS^J  + |T|2  E(S22)  T2Im  • (6.2-28) 

Consideration  is  nov  given  to  the  second  term, 

Re(EU*(k)  e2(k)  G'(k)}^  , (6.2-29) 

in  equation  (6.2-19).  Application  of  ccndition  five  in  section  6.1 
in  conjunction  with  the  development  in  appendix  D permits  the  preceeding 
term  to  he  expressed  as 

Re  ^E(e*(k)  e2(k)>  E(G'(k)}^  . (6.2-30) 

Equations  (6.2-1)  and  (6.2-2)  may  be  used  to  obtain  an  expression 

* 

for  e1(k)e2(k). 


e1(k)  e2(k)  « 


Jlm(T)  (xx 


(k)X2(k)  + Y1(k)Y2(k)j 


♦ Re(T)  ^(kjYgCk)  - X2(k)Y1(k)jJ 

|r«(T)  |xi(k)x2(k)  + Y^kjY^k)) 

- Im(T)  ^X1(k)Yg(k)  - X2(k)Y1(k))j 
E{e*(k)e2(k)}  = jlm(T)  ^(X^kJX^k) } + Ett^kjY^k) } j 

+ Re(T)  |E(X1(k)Y2(k)}  - Ett^kjY^k) 

|r  e(  T ) ^E(X1(k)X2(k)}  + EtY-^kjY^k)  }J 


- Irn(T)  ^E(X1(k)Y2(k)}  - E{X2(k)Y 


The  cross  correlation  terms  in  the  preceeding  equation  axe  evaluated 
in  appendix  E. 

Now  attention  must  be  given  to  the  evaluation  of  E{G'(k)}.  From 
equation  (6.2-3), 


G'(k) 


t-k-M 

Mk 


e^t) 

e2(l) 


SX1(« ) + JYX( 

±JT  (X2(*)  + i 

Uk-M  L 


(6.2-3** ) 
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k+M 

E 

t"k-M  I 
Mk 


[x1U)y2(D 


X2(£)Y1(Jl)]Be(T) 


[X1U)X2(£)  + Y1(£)Y2(£)]lin(T) 


|T | 2 [x|(£)  ♦ Y2U)] 


+J 


k+M 


lf<k 


[X1(t)X2(i)  + Y1(i)Y2(£)jRe(T)  - (X^jYgU)  - Xg(£  fcjU ) ]lm(T) 


|t|‘ 


[Xg(t)  + Yg( £) ] 


(6.2-35) 


The  assumed  uniformity  in  space  and  time  of  the  statistical  distributions 
allows  the  the  expectation  of  G'(k)  to  be  expressed  in  the  following 
form. 


E{G'(k) } 


(X1(£)Y2(£)  - X2(£)Y1(£)]Re(T)  + 

[X1( £ )X2( £ )+Y1(£  )Y2(  1 ) ]a»(T) 

[x|(£)  + y|(£)] 

E 


[X1(£)X2U)+Y1(£)Y2(£)]Re(T)  - [X^  £ )Y2(  O-XgUjY^  l)  ]lm(T) 

[x|(£)  + y2(£)] 


(6.2-36) 

Evaluation  of  E{G'(k)}  requires  that  the  expectation  of  the  quantities 

X1(£)Yg(£)  X2(f)Y1(£)  X1(£)Xg(£)  Y;l(£)Y2(£) 

X2(£)  + Yg(£)  * X2(£)+Y2(£)  ’ X2(£)+Y2(£)  * x|(0+y|(£) 

be  determined.  The  required  expectations  are  in  the  form  of  three  dimen- 
sional integrals , and  appendix  G shows  how  the  integrals  may  be  reduced 


to  two  dimensional  forms. 


The  third  term  of  equation  (6.2-19)  may  new  be  expressed  as 


E{|o'(k)|2}  * E{  Je-Ot)  J2> 


(2M+1) 


as  a result  of  the  work  culminating  in  equation  (6.2-30).  E{je2(k)|2} 

is  evaluated  in  the  same  manner  as  was  E(|e  (k)l2). 


E{|e2(k)|2}  « [E(X?(k)>  + E(Y2(k)}] 


(6.2-37) 


From  appendix  A, 


E{X2(k)}  - E{R^}  Ej^|*Sij  E{a2}  , and 

E{Y2(k) } - E{Hk) 


E{a2} 


(6.2-38) 


(6.2-39) 


Equation  (6.2-5)  may  be  used  to  show  that 


M2  * S21  + ,T|2  S22  +2Ia(TS21S22)  . 


(6.2-Uo) 


Thus, 

E(|a2|2}  * E(S21>  + |t| 2 E(S22)  + 2Io(T)  EtS^S^}  . (6.2-Ul) 

|G"(k)j2  is  indicated  in  equation  (6.2-17)  as  being  the  sum  of  two 
squared  quantities.  Each  of  the  two  squared  quantities  may,  through 
application  of  equation  (6.2-18),  be  shown  to  be  as  indicated  in  the 
following. 


*+M 


l-k-M 

Hk 


B«(T)^X2(A)Y1(i)~X1(t)Yg(Mj  + In(T)(kXgfit)Y1(Jt)-X1(i)Yg(t)) 


|T|2  (x|(*)  + Y2(*)j 


k+M  k+M 


E E 


Re(T)  (xg(t)Y1(t)~X1(&)Y2(X))  ■Hm(T)(xg(i)Y1(i)-X1(ft)Yg()t)) 


*»k-M  J«k-M 
A*k  J^k 

Hi 


|T|2  ^(1)  + Y2(*)j 


R«(T)  (xg(  1 )Y1(j)-X1(j)Yg(j))  + Im(T)(xg( J )Y1( 1 )-X^(j  )Yg(j  )) 
|T|2  ^(j)  + Y2(j)j 


(6. 2-U2) 


The  above  represents  the  first  term  of  equation  (6.2-17),  and  the  fol- 
lowing is  another  form  for  the  second  term  in  the  same  equation. 


k+M 


i«k-M 

i*k 


Re(T)  (xiU)X2U)+Y1U)Yg(l))+  3a(T)(x2(t)Y1(£)-X1U)Y2U)) 

!(*)) 


|t|2  ^x2U)  + Y2(i)'l 


k+M  k+M 


E E 


fc*k-M  J«k-M 

h*  jyk 
Hi 


Re(T)^X1(i)Xg(t)+Y1(i)Yj&)j+Im(T^X2U)^(Q-X1(fc)^uj) 


|T|2  ^X2(i)  + Y2(£ 


}J 


~R6(T)(xi(J)X2(3)+Y1(J)Y2(J))  + Xm(T)(xg(j)Y1(j)-X1(j)Y2(j)) 


|T|?  ^(j)  + Y2(,j)| 


(6.2-43) 


' 2 

Note  that  when  T is  real  the  following  simplified  form  results  for  |G'(k)|  . 


k+M 


i°'(k)i2  - y 


i=k-M 

Hk 


x2U)y1U)  - x1(h)y2(z) 


U)  + y2U)| 


2 


in  the  discussion  of  E{G'(k)}.  Evaluation  of  the  first  two  terms  of 
equation  (6. 2-1*5)  requires  knowledge  of 


(xg(i)Y1(t))  2 

^(OYgU))2 

T 2[X^(Jl)  ♦ Y2(£)]2 

• * 

T2[X2(£)  + Y2(£)]2 

E 


99 


| (xjUJXgU))2 

1 

, and  E 

(y*)Y2(.))2 

T2 

2;U)  + y2(a) 

2 

T2 

X§(£)  + Y22(l) 

2 

A technique  is  presented  in  appendix  F whereby  the  above  expectations 
may  be  reduced  to  two  dimensional  integrals. 


6.3  Lover  Bound  on  Performance 

Conditions  were  specified  in  chapter  U under  which  the  lower  bound 
of  performance  resulted  for  the  ordinary  system.  Those  same  conditions 
will  be  applied  here  to  the  cancellation  system  equations  for  the  purpose 
of  comparing  the  performance  of  the  cancellation  system  with  the  known 
lower  bound  of  performance  of  the  ordinary  system.  The  conditions  de- 
scribed in  chapter  U require  that 


2 I “ 1 

cos  (</')  =l  1 . 

\ r J 

(6.3-1) 

/v  \2 

? / Voo  1 

sin  (ijz)  ; 1 - y-  » 

\ rl 

(6.3-2) 

sin(c)  * 1 > and 

(6.3-3) 

cos(c)  * 0 

(6.3-M 

Under  these  conditons,  the  S matrix  elements  become 


hi  - (gx-gy)  (V-/V2 

S12  * (6,-gy)  <V„/Vr>  Vl-(v„/vr)a 


s,2  --  V«y>  (V-/Vr>  * *1 


(6.3-5) 


(6.3-6) 


(6.3-7) 


fa."1 


Other  results  from  chapter  four  which  will  be  used  here  are 
2 


2 

a » 
®x 


[jjts]  d°  f8'25  ' 0-2U1,do  * °-023doJ  ’ 


( 6. 3—8) 


2 

®y 


L-rrf  < '*■ 


33  + 0. 112df ] , 


(6.3-9) 


■ trfel 


d®  [0.0024  + 0.00l*15d2] 


(6.3-10) 


gx*/  "[ex  10^]  ^ “ °,0672d0  + 0*0115do  " 0*002075dJ]  , (6.3-11) 


-Vs 

V,  = 10.105  [1-e  0 ] , 


(6.3-12) 


and  the  drop  size  distribution  is 


P(dQ)  » X-  d2  e Y 0 


(6.3-13) 


Evaluation  of  the  cancellation  system's  output  equation  requires 
v \owledge  of  various  moments  of  the  S matrix  elements.  The  required 
moments  are  given  below. 


E{Sn}  * E{«x}  - 2E{exVgy)(VV2}  + E{(gx“g  )?(v«/v  >**>  (6.3- 


y * r 


H*) 


E(S22)  » E((gx-gy)2  (Vw/Vr)2}  - E{(gx-g )2  (VM/Vr)U) 


(6.3-15) 


’TPPlSPHW-.'J.'iJ1  unimiim  ■ i.»iipn. «,',■■•'  >.i»f iipm.*..' i . t - 
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E{sl2)  * E{(gx-gy)2  (VjVr)k)  + 2E{gy(gx-gy)  (V^)2}  + E{g2}  (6.3-16) 


E{S11S12}  * E{8x(gx-Sy)  (V«/Vr))  “ <|>  ^(g^)  (V^)3} 


- E{(gx~gy)2  (V^/Vj.)3}  + (|)  E{(gx-gy)2  (VM/Vr)5}  (6.3-17) 


E{S11S22}  " E{gx(gx“gy)  ^VVr)2>  - ^(«x-gy)2  (V„/vr)S 


+ E{gxgy)}  - E{gy(gx-gjr)  (Va/Vr)2) 


(6.3-18) 


Evaluation  of  equations  (6.3-11*)  through  (6.3-18)  is  accomplished 
in  appendix  G,  and  the  following  material  shows  how  to  apply  the  results 
obtained  to  determine  the  output  noise  for  the  cancellation  system. 

The  following  assumes  that  T is  real. 

[j  p. 

E{  |e^(k)  |2}  = |A|2  E(Nk)  E{°  • ps^}  + T2  EtS^jJ  (6.3-19) 


E(e1(k)e2(k) } = 


(E{V 


k)Y^(k)}  - E{X2(k)Y1(k)>^ 


-J 


^E{X1(k)X2(k)}  + E(Y1(k)Y2(k):j 


E +i|LTE{Nk}  E{GU(a,0)}  pTjE^S^-S^} 


(6.3-20) 


-j  E(SllS?1  * S12S2/J] 


(6.3-21) 


J|  — •» 

E{  |e2(k)|2}  * 72  |A|2  E(Nk)  E{°- j,-^?  * E{ ) + T;>K(Sp?)J  (6.3-22) 


rain  intensity  (mm. /hr. 
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Fig-  6.3-5. 

Lower  Bound  on  Cancellation  Improvement 
Cancellation  System  vs.  Ordinary  C.P.Syst 


Fig.  f.3-6. 

Lower  Bound  on  Cancellation  Improvement 
Cancellation  System  vs.  Ordinary  C.P. System 
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The  terms  presented  in  equations  (6.3-19),  (6.3-21),  and  (6.3-22) 
may  be  reduced  to  simple  algebraic  expression  by  the  techniques  presented 
in  appendix  G,  but  no  means  has  as  yet  been  found  that  accomplishes  the 
same  reduction  for  the  moments  of  G'(k).  Thus , the  moments  of  G ' ( k ) 
have  been  evaluated  by  numerical  integration  techniques  presented  in 
appendix  H. 

When  the  numerical  values  of  equations  (6.3-19),  (6.3-21),  and 
(6.3-22)  are  combined  with  the  evaluated  moments  of  G'(k)  according 
to  equation  (6.2-19)  the  result  is  an  expression  for  the  output  noise 
power  of  the  cancellation  system.  The  value  for  noise  power  thus  obtained 
may  be  taken  in  ratio,  according  to  equation  (5.^-5),  with  the  output 
power  of  the  reference  system  to  obtain  the  cancellation  ratio.  Note 
that  the  reference  system  output  power  is 

E{  |er(k)REF|2}  = |A|2  E(Nk)  E{GU(a,B))  EtS^}  . (6.3-23) 

Numerical  data  are  presented  in  the  following  figures  which  describe 
the  performance  of  the  cancellation  system  relative  to  the  reference  system 
for  various  rain  sources  over  a wide  range  of  rain  intensity. 

6,h  Upper  Bound  on  Performance 

As  indicated  in  chapter  four,  the  conditions  which  lead  to  the  upper 
bound  of  performance  for  the  ordinary  system  are 

sin‘  (0  = 0 , ( 6.  h-1  ) 

eos2(0  = 1 , 


( C.l-2 ) 


The  S matrix  elements  are  defined  to  be 


SXl'<V 


(6. It- 


S^2  ■ 0 , and 

V 

S22  * «y  + (g*V  ^ 


(6.  U— 
(6.U- 


under  these  constraints.  Evaluation  of  the  cancellation  system  output 
equation  is  carried  out  in  the  same  manner  as  was  used  in  the  analysis 
of  the  lower  bound  case  with  the  exception  that  the  above  expressions 
for  the  S matrix  elements  are  used  rather  than  those  given  in  section 
6.3. 

Numerical  data  in  the  form  of  the  cancellation  ratio  is  shown  in 
the  following  figures.  Again  various  rain  sources  are  considered,  and 
the  variation  in  performance  with  rain  intensity  is  also  shown. 


Upper  bound  of  cancellation. 


Fig.  6.^-2.  Upper  Bound  on  Cancellation  Improvement 

Cancellation  System  vs.  Ordinary  C.P.Syatem 
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7.°  CONCLUSIONS 

Based  on  the  computer  simulations  and  the  rain  model  described 
in  paragraph  6*1  of  this  report,  rain  cancellation  improvement  of 
1 dB  to  over  8 dB  can  be  achieved  over  that  of  circular  polarization 
when  a rain  polarization  estimater  is  used  to  determine  the  proper 
null  polarization.  These  results  are  illustrated  in  Figures  6.3-4 
through  6.3-6  and  6.4-2  and  represent  rain  of  certain  origins. 

A study  of  this  tvpe  must  be  followed  by  an  experimental  program 
in  order  to  subject  the  technique  to  real  world  environments. 
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APPENDIX  A:  ANALYSIS  OP  MOMENTS  FOR  RANDOM  SUMS 


Asssume  that  it  is  desired  to  find  the  statistical  moments  of  the 


random  quantity 


X - [ Uj 


where  n and  are  "both  random  in  nature.  Beckmann  [46]  has  shown  that 
the  density  function  for  the  variable  X i3 


f(x)  * l P(n)  p(X'n  ) 


where  P(n)  is  the  discrete  density  function  of  n and  p(XVn;)  is  the  con- 
ditional density  function  of  X given  a specific  value  of  n. 

The  mean  value  of  X is  calculated  as  follows. 


E{X)  « x f(x)  dx 


oo  00 

| x l P(n 


) P(x'n/)  dx 


00  oo 

1 P(n)  | X!n)  p!X(n))  dx(n> 


n«0  -« 


is  Gaussian  then 


’,-"T  ■•"ni"  -.tww 


random.  Beckmann  [U6]  observes  that  terms  having  the  form  exhibited  by 

y1(k),  Xg(k),  and  y^(k)  tend  to  be  Gaussian  in  nature  if  the  number 
of  terms  in  the  appropriate  summations  is  large.  The  number  of  terms  in 
the  summations  corresponds  to  the  number  of  rain  particles  in  a range  cell, 
and  the  material  of  Chapter  3 indicates  that  the  number  of  particles  in  a 
range  cell  is  generally  very  large.  Thus,  the  density  function  for  x^(k)  is 


p(x  (n,(k))  = -i-  exp  i 

cr/sFirn 


-(x.  ^(k)  - na)^ 


(A-16) 


where 


Eju2(ait  Si)“11(k)  cos(Yli(k)  - 0^ 
E{G"^ai*  6i^ali^  cos(Yli(k))cos(6i)^ 
Ejo2(ai»  8i)c‘11(k)  sin(Y11(k))sin(6i)| 


(A-17) 


(A-18) 


o2  * E{(U4  - a)2} 


(A-19) 


The  condition  of  uniform  particle  distribution  means  that  6^,  as 
defined  in  equation  (5.2-6),  is  uniformly  distributed  between  the  symmetric 
extremes  + 0,.  and  - where 

M M 


0M  “ 260AZM 


(A-2C) 


Thus,  the  density  function  for  0^  is 


s<V  ■ TT 


(A-21) 


w 


If 


t 


\ 

f 


9M=±mir  , m = 1,  2,  3,  ... 


( A-22 ) 


then  a in  equation  (A-1T)  equals  zero,  but  another  statement  of  the  conditior 
given  in  equation  (A-22)  is 


ct  9 

2 Azm  = -f = Bj 


( A-23 ) 


where  tp  is  the  transmitter  pulse  length  and  o is  the  speed  of  light.  Equation 
(A-23)  may  be  rearranged  to  obtain 

t = ±m/fQ  (A-2U) 

where  fQ  is  the  carrier  frequency  of  the  transmitter.  In  the  following 
work  it  will  be  assumed  that  the  condition  of  equation  (A-2t)  is  satisfied. 

Application  of  conditions  one  and  two  as  stated  in  section  6.1  to 
equation  (A-l8)  and  (A-19)  gives  the  following. 


2 

a = E(G  (a,  6)0^  cos  y}  E{cos  9} 

+ E(G2(a,  6)01^  sin  y)  E{sin  9}  . (A-2 S) 

= 0 (A-2 6) 

a2  * E(Ui2}  (A-2T) 

= E{GU(a,  6)0^  cos2(y1  - 8)}  (A-20) 

= (~)[E(GU(a,  B)^2} 

+ E{G1'(a,  8)a12  ccs(2y1  - 29)}]  (A-fo) 

» (“)E{G1*(a,  P)>  E{ a}2 } (A-3C) 
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Thus,  from  equations  (A-10)  and  (A-15) 


E{x1(k)>  * 0 

E{x12(k)}  * E{Nk)  E{-G-^»  E{ai2} 


(A-31) 


(A-32) 


The  moments  for  Xg(k),  y^k),  and  yg(k)  are  evaluated  in  a manner 
similar  to  that  used  above  for  x^(k),  and  the  results  are  presented  below. 


E(x2(k)}  * 0 


(A-33) 


E{x22(k)}  « E{Nk>  e{o22> 


Ety-^k)}  = 0 


(A-3U) 


(A-35) 


E{yi2(k)}  * E{Nk}  E{P--^g--6^)  E{ai2} 


(A— 36) 


E{y2(k)}  a o 


(A- 37) 


rG^(o,  8) 


E{y0  (k) } a E{N.  } E{"  E{a  } 


(A-38) 


122 


APPENDIX  B:  THE  MEAN  VALUE  OF  N, 

k 

The  material  in  Chapter  3 indicates  that  the  volumetric  particle 
density  is  of  the  form 

K? 

PD  * K I (B-l; 

where  I is  the  rain  intensity,  and  and  are  constants  determined 
"by  the  type  of  rain.  Additionally,  1 may  he  represented  as 

I « I + 61  (B-2) 

where  1^  is  the  mean  value  of  the  rain  intensity,  and  61  is  the  random 
component  of  the  rain  intensity. 

Calculation  of  the  total  number  of  particles  is  accomplished  by 
forming  a product  of  the  cell  volume,  V , and  the  particle  density,  PD. 

it 

N = PD  • Vk  (B-3) 

( B-L ) 

(B-5) 

is  a zero  mean  Gaussian  variate  and  that  most 

(B 

, and  ( B-T ) 

(H-fl) 


= hU0  * 61)  Vk 


K,  (I  K, 

’Wo 2 2 


It  will  be  assumed  that  61 
of  the  time 


(fi>  « i. 
0 


Thus 


Kp  K 61 

*k  = Wo  [1  + \ ] 


K 1 N } = K.,  V I 


K, 


1 h 0 


r1' 


The  values  of  and  Kg  are  provided  in  Table  B.l. 


Table  B.l.  Parameters  K and  K for 
severed  rain  sources . 


SOURCE  OF 
RAIN 

K1 

K2 

MELTING 

HAIL 

381+ 

0.31 

MELTING 

SLEET 

203 

0.31 

MELTING 

SNOW 


87-5 


0.39 
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APPENDIX  C : EVALUATION  OF  THE  FOURTH  MOMENT  OF  THE  RADIATION 

CHARACTERISTIC 

The  radiation  characteristic  assumed  for  this  work  is  in 


any  plane  passing  through  the  axis  of  propagation 


G(x.n)  « 


2J2(x) 


= 3.83  (18) 


where  Jv(x)  is  the  Bessel  function  of  the  first  kind  and  first 
order,  and  where  x and  ft  are  defined  in  FIG.  C-i.  It  is  assumed 
that  the  particle  position  parameters,  x and  0,  are  statistically 
independent  of  one  smother  and  that  the  distributions  of  these 
variables  are  uniform. 


The  fourth  moment  of  G(x,ft)  which  is  required  for  the  analysis 
in  chapter  5 is  of  the  form 


X 2tt 
*u 


‘(X.s»)  - f } PV»r 

0 0 -1—  x -- 


dfi  dx 


!6jJ(x) 
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V X ■* 


16  [V<‘> 
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Attempts  at  evaluating  equation  (C-5)  analytically  in  closed 
form  have  thus  far  been  unproductive,  and  therefore,  an  evaluation 
was  affected  by  means  of  numerical  integration.  The  numerical 
method  used  was  a Newton-Cotes  five  point  formula  as  suggested  by 
Weeg  and  Reed  [1*8]  with  an  increment  in  the  independent  variable, 
X,  of  0.01. 

E {G1*(X,0)}  ■ 0.70808 


(C-7 ) 
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APPENDIX  D:  STATISTICAL  CONSIDERATIONS  OF  A PRODUCT 

It  is  assumed  that  several  random  processes  exist  and  may  be  divided 
into  two  groups,  x and  y,  such  that  the  members  of  the  groups  exhibit 
the  characteristics  that  x^  and  x^  are  statistically  independent  if  i ^ J, 
that  y^  and  yj  are  statistically  independent  if  i ^ J , and  that  x..  and  y^ 
are  statistically  independent  if  i ^ J.  Under  the  stated  conditions,  one 
may  represent  the  Joint  density  function  of  elements  of  the  two  sets  as 
follows  by  application  of  the  concept  of  conditional  probability. 

p(x. , . . . , x , y , . . . y ) = 
x n i n 

p(x1,y1)p(x2,  ...  xn,  y2,  ...  yn  |xlt  y^)  (D-l) 

The  assumed  independence  requires  that  the  individual  or  Joint  out- 
comes of  x1  and  y.^  in  no  way  influence  the  realizations  observed  for  the 
other  x and  y elements.  Thus 


p(*2»  .. 

*•  Xn’ 

y2.  • • • 

* yn 

I*!* 

h1  = 

p(*2»  * * 

xn. 

y2,  • • • 

yn}’ 

and 

(D-2) 

p(xx»  •• 

xa. 

yx»  • • • 

• yn) 

= 

p(xx.  y± 

)p(x2» 

. . . , X 

n 

t y2, 

• ' • » 

»«>■ 

(D-3) 

Let 

a function  f 

exist 

such 

that 

f * x^  • 

yl  ‘ 

g where 

CD— ) 

8 * g(*2 

• • • • » 

V y2» 

• » • 1 

ynK 

(D-5) 

Then 

Elf)  * / ...  /x1y1g(x2,  ....  yn)p(xlf  x^). 
2& 

p(x0,  ....  y )dx  ...  dy 
d n l n 


(D-6) 
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The  following  expressions,  cttained  from  equations  (6.2~k)  and  (6.2-5) 
may  be  applied  to  equation  (E-12), 


cosy.. 


Su  ± a12  Irn(T) 


(E-13 


siny1 


±s12  Re(T) 


(E-l 


cosy. 


S21  4 S22 
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Thus  equation  (E-12)  may  be  rewritten  as 


E{xx(k)x2(k)}  = E(Nk)  E(G  (o,B)J 
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E(~ — coby1cosy2 
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The  fourth  term  to  he  evaluated  is  E(Y^(k)Y2(k)}. 
Y1(k)Y2(k)  = 3i)ali(k)sin^li(k)-6ij 


c2(“i>e1)“2iCk’sln  (f2i(lt)-e1) 


(£-19) 
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APPENDIX  F:  THE  REDUCTION  IN  DIMENSION  OF  INTEGRALS  USED  TO  EVALUATE 

THE  EXPECTATION  OF  G'(k)  AND  |G'(k)!2 

Evaluation  of  E{G'(k)}  requires  the  determination  of 


Jx1(t)x2(t) 

|x^U)  + Y2U) 


, » 


Y1(£)Y2(£) 
X2(£)  + y|(£) 


, E< 


^(OYgU) 

|x^(£)  + y|(£) 


, and  E< 


|X2(£)Y1(£) 

| X^U)  + Y2U) 


E< 


Ix^OXgU) 

|X2U)  + T2U) 


(i) 

Y2(£) 


p(xi(»),X2(t),Y1(£),Y2UjdX1(£)...dY2(A) 


(F— 1 ) 


Wilks  [20]  indicates  that  integration  of  the  density  function,  p(*), 
with  respect  to  Y^A)  will  result  in  a marginal  density  function  of 
three  variables  with  the  same  means  and  variances.  Thus, 


K(£)X2U) 


(x|(£)+Y2(£) 


fffxAi)x(i)  , x 

“I  I / o H P (X,(£),X  (£),Y0(£faX  (£)... dY.(£) 

JJJ  X2(£)+Y2(£)  V 1 ’ 2 " / 1 2 


(F-2) 


Appendix  A points  out  that  the  X and  Y parameters  are  zero  mean 
Gaussian  variates,  and  thus  the  Joint  density  function  may  be  expressed 
in  the  form 


/ x -[c,x2  + C X + c3] 

p^UhXgUJ.Y^JDJ  = l^e  1 1 21 


(F-3) 


The  quadratic  form  appearing  in  the  exponential  is  known  to  be  positive 
definite  [1*7],  and  thus  C^  is  known  to  be  positive  [50].  Therefore, 
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3"  5c^J 


dX2U)dY2U) 


(F-U ) 


where  the  inverse  of  the  covariance  matrix  of  the  variables  X^(H), 
Xg(fc),  and  Y2(£)  is 


a b c 
d e f 
jg  h i 


(F-5) 


and 

C;L  = a/2  (F-6) 


c2  * [(b+d)X2(t)+(c+g)Y2U)]/2 


(F-7) 


c3  = [eX2U)+(f+h)X2U)Y2U)+iY2(£)]/2  (F-8) 

K = -7-i — . (F-9) 

^ (2w)3/2  • ^Tet (y 

A similar  approach  is  used  in  evaluating  the  other  expectation 
terms.  The  results  are  presented  below. 


where  the  inverse  of  the  covariance  matrix  of  the  variables  Y^fc), 
Y^U),  and  X2U)  is 


PfpBWJFpfW  ■.  mu  W '.g  jf  t vwi  wjij  piwjm  mm  iniUcpii 

ar 
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a;1- 


a b c 
d e f 
[g  h i 


(F-ll) 


and 


C.  * a/2 


CK  * [(h+d)Y0(£)  + (C+g)X0(£)]/2 


C6  = [eY2(£)  + (f+h)X2U)Y2U)  + ix|(£)]/2 


Kn  = 


2 ,372 


(2ir) • /det(AF) 


4-oo 


(F-12) 


(F-13) 


(F-lU) 


( F-l 5 ) 


E< 


|X1(A)Y2(£) 

|X2(£)+Y2(£) 


-K,C0/tT\  /yo(£)  \ -[Crl~-) 

• e 1 


12 


20^1  I \x2(£)+y|(£); 


dX2(£)dY2(£) 


(F-16) 


+00 


E< 


|X2(£)Y1(£) 

|X^(£)+Y2(£) 


'-KgC^/jT^  ^X2(£) 


2c2^2  /\x2(£)+y?u) 


-[Cr-  rl— ] 

j 6 4CU  dX2(£)dY2(£) 


( F-1T ) 


The  following  expressions  are  derived  for  use  in  evaluating 
E{  |0'(k)|2>. 


( £ )X_ ( £ ))  ^ /*/■/■  X?(Jl)x|(£)  \ 

\ -///  ~~2 ^7772?  X1(0,X  (£).Y 2(£)'dX  (£)... dY  (£) 

*U)+Y2(£)]2  JJ7  [X^(£)<(£)]2  V 1 2 2/1 


|[x|(-,.*g 


(F-l8) 
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APPENDIX  G:  MOMENTS  OF  THE  S MATRIX  ELEMENTS 

Expectations  of  Sll*  S12*  S22*  SilS12’  S11S22*  and  are 

required  for  the  evaluation  of  certain  terms  appearing  in  ehpater 
six.  These  terms  will  be  evaluated  in  the  stated  order  here. 


8.25  x 81  _ 0. 2kh  x 101  + 0.023  * 111' 
(•|+y)9  (§tY)n  (-^ty)12 


8.25  x 8|  _ 0.2U*  x ioi  + 0.023  x 131 
(l+y)9  (l+y)11  (l+y)12 


fWBWWWIfWWWi  gw  mn.iuw1  wiw  ipmm  '»■<»■ 


rrn-*?f>r.r - 
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1 

re. 29  x 8! 

0.0672  x 10! 

. 0.0115  x 111 

0.002075  x 12! 1 

i 

[ T* 

Yn 

y12 

Y13  J 

o 

["8.29  x 8! 

0.0672  x 10! 

, 0.0115  x 11! 

0.002075  x 12 jl 

c. 

L«h>9 

(“*Y)U 

(|+Y)12 

(^>13  J 

Ts 

.29  x 8i 

0.0672  x 10!  , 

0.0115  x U! 

0.002075  x 12! 

L (i+y)9 

(i+y)11 

(i+y)12 

(i+y)13  [ 

(G-6) 


.1 

j 0.0024  x 10! 

, 0.00415  x 12! " 

, To. 0024  x 10! 

| 0.00415  x 12! 

i 

L Y11 

y 1 3 

L (lfY)11 

(|+y)13 

[*0.0024  x 10! 

( 0.00415  x 12! 

. [*0.0024  x 10! 

( 0.00415  x 12! 

▼ 0 

[ ( 1+Y ) 1 1 

( 1+Y ) 1 3 

(fv>13 

♦ 

0.0024  x 10! 

> 

, 0.00415  x 12 1*1 

1 

(0-7) 

. 

(2+Y)11 

(2+y)13  j 

f 
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The  following  notation  is  adopted  in  order  to  simplify  represen- 
tation of  the  following  work. 


\ = Y3/2 

A12  = 0.33x8!, 

G1C  = (1+Y)11 

a2  = 8.25x81 

A13  = 0.112x10! 

G^  = (1+y)12 

a3  = 0.21414*101 

G1  = Y9 

°12  = (1+Y)13 

\ = 0.023X11! 

g2  = y11 

°13  * 

A5  = (I0.105/Vr)2 

g3  - y'2 

6u,- 

A6  = 8.29x8! 

% = Y13 

°i5  * {¥y)X7 

A^  * 0.0672x101 

3,  - (~YY>9 

C16  ■ <^>13 

Ag  = 0.0115x11! 

°6  ■ <h>U 

G1T  a (2+y)9 

A?  * 0.002075x12! 

°7  = (2*t)12 

Gi8  a ( 2+Y ) 1 1 

A^  » 0.002UxiO! 

c8=  (|n)>3 

G19  = (2+y)12 

AX1  = 0.00)415x12! 

Gg  = (1+y)9 

G20  a (2+y)13 

°21  " ^Y)11 

°?2  - 
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fv 

\2  | » 


3!2  ■ (^V‘  Ivfl  - <*  V ft] 


(G-8) 


E(S^>  - E^-g/ 


E<(gx“S}  w 


'll 


10  12J 


2 hi 

E^g^-Sy)  (Vw/Vr)  | iG  evaluated  in  equation  (G-7). 


\2  V« 


22 


(«x-«y)  \yrj  + V^’V^rj  +gy 


(G-9) 


(G-10) 


(G-ll) 


E{S22>  » E^-gJ2 


2E  <&. 


r(«bc-ey)[^r 


) + E{g^} 


(G-12) 


'^2  . hi 

°13  °U 


(G-13) 


y (g-iu) 


( G—  22 ) 


(G-23) 


(G-2U) 


(G-25) 
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APPENDIX  H: 


NUMERICAL  INTEGRATION  TECHNIQUES 


Numerical  integration  will  be  used  to  evaluate  the  required  moments 
of  G'(k).  The  work  of  Appendix  G has  shown  that  the  desired  moments  may 
be  expressed  as  two  dimensional  integrals  extending  over  (-«,  +») , but 
for  the  purposes  of  numerical  integration  it  is  desirable  that  the  range  of 
integration  by  {0,  l).  The  change  of  integration  range  is  accomplished  in 
the  following  manner. 


i 


-H»  0 00 

| f(x)  dx  * [ f(x)  dx  + | f (x)  dx 

-CO  _oo  0 

In  the  range  (-«,  0),  let  x = In  y so  that 

0 1 

f ( x ) dx  = f f(in  y)  ^ , 

. j y 

-00  0 


(H-l) 


(H-2) 


and  in  the  range  (0,  +°°),  let  x = -In  y'  so  that 


0 1 

f(x)  dx  = - [ f(-*n  y)  = f(-*n  y') 

j y y 


Thus , 


' r(x)  ax  - j tf<*°  ?>  t d=ia  ar)l  dy. 


( H— 3 ) 


(H-U) 


This  type  of  transformation  is  suggested  by  Davis  and  Rabinowitz  [51]. 

The  actual  numerical  integration  scheme  is  basically  a Monte  Carlo 
technique,  but  instead  of  using  a multidimensional  random  variable  generator 


to  provide  the  independent  variables,  a set  of  linearly  independent 
equidistributed  sequences  is  used.  A justification  for  and  iiscussion 
of  this  technique  is  provided  by  Davis  and  Rabinowitz  [5i].  Integration 
was  carried  out  using  one  thousand  points  of  evaluation  in  each  quadrant 
of  the  independent  variable  space,  and  checks  on  the  accuracy  of  integration 
were  developed  as  indicated  in  the  following  figures. 


! 


Fig.  H_2.  Accuracy  of  numerical  integration. 
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PAX  i»  the  principal  AFX  organization  charged  with 
planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  information  sciences.  Intelli- 
gence, command,  control  and  communications  technology, 
product*  and  services  oriented  to  the  needs  of  the  USAF. 
Primary  XADC  mission  areas  are  communications , electro- 
magnetic guidance  and  control,  surveillance  of  ground 
and  aerospace  objects.  Intelligence  data  collection  and 
handling,  information  system  technology,  and  electronic 
reliability,  maintainability  and  compatibility . PAX 
ban  mission  responsibility  as  assigned  by  AFX  for  de- 
monstration and  acquisition  of  selected  subsystems  and 
systems  in  the  intelligence , stepping,  charting,  command , 
control  end  communications  areas. 


